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ABSTRACT 

We analyze the absorption and emission-line profiles produced by a set of simple, cool gas wind 
models motivated by galactic-scale outflow observations. We implement monte carlo radiative transfer 
techniques that track the propagation of scattered and fluorescent photons to generate ID spectra 
and 2D spectral images. We focus on the Mg II AA2796, 2803 doublet and Fe II UVl multiplet at 
A w 2600A, but the results are applicable to other transitions that trace outflows (e.g. Na I, H I 
Lya, Si II). By design, the resonance transitions show blue-shifted absorption but one also predicts 
strong resonance and fine-structure line-emission at roughly the systemic velocity. This line-emission 
'fills-in' the absorption reducing the equivalent width by up to 50%, shift the absorption-lin centroid 
by tens of km s~^, and reduce the effective opacity near systemic. Analysis of cool gas outflows 
that ignores this line-emission may incorrectly infer that the gas is partially covered, measure a 
significantly lower peak optical depth, and/or conclude that gas at systemic velocity is absent (e.g. 
an interstellar or slowly infalling component). Because the Fe II lines are connected by optically- 
thin transitions to fine-structure levels, their profiles more closely reproduce the intrinsic opacity of 
the wind. Together these results naturally explain the absorption and emission-line characteristics 
observed for star- forming galaxies at z < 1. We also study a scenario promoted to describe the outfiows 
of z ~ 3 Lyman break galaxies and find prfiles inconsistent with the observations due to scattered 
photon emission. Although line-emission complicates the analysis of absorption-line profiles, the 
surface brightness profiles offer a unique means of assessing the morphology and size of galactic-scale 
winds. Furthermore, the kinematics and line-ratios offer powerful diagnostics of outflows, motivating 
deep, spatially-extended spectroscopic observations. 
Subject headings: galaxies: formation - galaxies: starbust 



1. INTRODUCTION 

Nearly all gaseous objects that shine are also observed 
to generate gaseous flows. This includes the jets of pro- 
tostars, the stellar winds of massive O and B stars, the 
gentle Solar wind of our Sun, the associated absorption 
of bright quasars, and the spectacular jets of radio-loud 
AGN. These gaseous outflows regulate the metal and 
dust content and distribution within the objects and 
their surroundings, moderate the accretion of new mate- 
rial, and inject energy and momentum into gas on large 
scales. Developing a comprehensive model for these flows 
is critical to understanding the evolution of the source 
and its impact on the surrounding environment. 

Starburst galaxies, whose luminosity is dominated by 
H II regions and massive stars, are also observed to drive 
gaseous outflows. These flows are generally expected 
(and sometimes observed) to have multiple phases, for 
example a hot and diffuse phase traced by X-ray emission 
toge ther with a cool, denser phase traced by Ha emission 



(e.g.l Heckman et a l.ill990HMartinlll999l: [Strickland et al.l 
|20M Kobulnicky fc Martini 120101 ). Several spectacular 



examples in the local universe demonstrate th at flows can 
extend to up to ^ 10 kpc fr o m the galaxv (.Lehnert et ahl 
119991 IVeilleux et all l2003t iWestmoauette et al.l \2Qm 
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carrying significant speed to escape from the gravita- 
tion al potential well of the gala xy's dark matter halo 
(e.g. [Strickland fc Heckmanll2009l ). 

Galactic outflows are also revealed by UV and optical 
absorption lines, e.g. Na I, Mg II, Si II and C IV transi- 
tions. With the galaxy as a backlight, one observes gas 
that is predominantly blue-shifted which indicates a flow 
toward Earth and away from the galaxy. These transi- 
tions are sensitive to the cool (Mg II; T ^ lO'^K) and 
warm (C IV; T ^ lO^K) phases of the flow. The inci- 
dence of cool gas outflows is nearly universal in vigor- 
ously star-forming galaxies; this includes systems at low 



z which exhibit Na I an d Mg II absorption (Rupke et al 
2005a; Martin 2005; Sa to et all [2001 [ Martin fc Bouche 



i2009i : iChen et al.,,2010. ) , z ~ 1 star-forming galaxies with 
winds traced by Fe II and Mg II transitions (jWeiner et al.l 
l2009HRubin et al.ll2010bD . and z > 2 Lyman break galax- 
ies (LBGs) that show blue-shifted Si II, C II, and O I 
transitions ( Stcidcl et all 119961: ILowenthal et all 119971: 
[Shaplev et al,,,2003. ). 

The observation of metal-line absorption is now a well- 
established means of identifying outflows. Furthermore, 
because the X-ray and Ha emission generated by winds 
is faint, absorption-line analyses have traditionally been 
the only way to probe outflows in distant galaxies. How- 
ever, very little research has been directed toward com- 
paring the o bservations agains t (even idealized) wind 
models (e.g. iFuiita et al.l |2009() . Instead, researchers 
have gleaned what limited information is afforded by 
direct analysis of the absorption lines. The data di- 
agnose the speed of the gas relative to the galaxy, yet 
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they poorly constrain the optical depth, covering frac- 
tion, density, temperature, and distance of the flow from 
the galaxy. In turn, constraints related to the mass, 
energetics, and momentum of the flow suffer from or- 
ders of magnitude uncertainty. Both the origin and 
impact of galactic-s cale winds, therefor e , remain open 
matt e rs of debate (iMurrav et al.l 120051: iSocrates et al.l 



TABLE 1 

Mg II AND Fe II Transitions Considered 



200a [Strickland fc HeckmanI l2009t iSimcoe et al.l 120021: 



Qppenheimer fc Davell2006t iKeres et al. 

Recent studies of z < 1 star-forming galaxies have re- 
vealed that the cool outflowing gas often exhibits signif- 
icant resonant-line emission (e.g. Mg II, Na I) in tan- 
dem with the nearly u biquitous blue-shifted absorption 
JWcincr ct al. 2009:' M artin k Bouchell2009HRubin et al.l 
[2010bi : iChen et al.n2010l ). The resultant spectra resemble 
the P-Cygni profile characteristic of stellar winds . This 
phenomenon was first reported by iPhillipsI ()1993f ) , who 
observed blue-shifted absorption and red-shifted emis- 
sion for the Na I transition in the spectru m of the local 
starbu rst galaxy NGC 1808. More recentlv. I Weiner et al.l 
(|2009t ). who studied Mg II absorption in z ~ 1.4 galaxies, 
reported Mg II emission in a small subset of the indi- 
vidual galaxy spectra of their large sample. These were 
excluded from the full analysis on concerns that the emis- 
sion was related to AGN activity. The stacked spectra 
of the remaining galaxies, however, also indicated Mg II 
emission, both directly and when the authors modeled 
and 'removed' the v ~ Qkm s~^ absorption component. 
The authors suggested the emission could be related to 
back-scattered light in the wind, but presumed that it 
was related to weak AGN activity. Ve r y simi lar Mg II 
emission was observed bv iRubin et al.l feoiObl) who re- 
peated the analysis of iWeiner et al.l ( 2009() on a set of 
lower redshift galaxies. 

Bright Mg II line emission has al so been reported for 
individual galaxies at z ^ 0.7 by (jRubin et al.l l20l"0aL 
see also Rubin et al 2011, in pr ep.). In the i r anal- 
ysis of a single galaxy spectrum, IRubin et al.l (|2010aD 
further demonstrated that the Mg II line emission is 
spatially extended, and used the size of the emission 
to infer that the wind extends to at least 7kpc from 
the galaxy. These authors additionally detected line 
emission from non-resonant Fe II* transitions, and at- 
tributed the emission to fluorescence powered by Fe II 
resonant absorption. In other words, these photons are 
re-emitted by the wind into our sightline, and are analo- 
gous to the emitted photons in a P-Cygni profile. Line- 
emission that may be related to outflows is also ob- 
served for z ~ 3 LBGs in the reson ant Lyg transition 
and non- resonant S i II* transition^ (jPettini et al.l 120021: 
IShaplev et ani200l . This emission likely arises from a 
completely different physical process than those generat- 
ing X-ray and Ha emission (e.g., shocks), and presum- 
ably probes both the cool gas at the base of the wind and 
the outskirts of the flow (i.e., wherever a given transition 
is optically thick) . A comprehensive analysis of the scat- 
tered and fluorescent emission related to galactic-scale 
outflows ( e.g. via deep integral-field-unit [IFU] ob serva- 
tions, Sw inbank et al.ll2005l : IWestmoquette et al.ll2008| ) 
may offer unique diagnostics on the spatial extent, mor- 
phology, and density of the outflow from distant galaxies, 
eventually setting tighter constraints on the energetics of 

* See also [France et aLl l|2010l ) for a new example at z ~ 0. 
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the flow. 

Although astronomers are rapidly producing a wealth 
of observational datasets on galactic-scale winds, a key 
ingredient to a proper analysis is absent. Just as compar- 
isons between observed supernova lightcurves and spec- 
tra and radiative transfer calculations of of SNe models 
have provided crucial insight into the physics driving, 
e.g., standard candle relationships for both Type la and 
IIP SNe ( K asen fc Woos lcv 2007, 2009), modeling of the 
observable signatures of galactic outflows is necessary for 
understanding the physical properties of the gas and ul- 
timately the physics driving the flows. In this paper, 
we take the first steps toward modeling the absorption 
and emission properties of cool gas outflows as observed 
in one dimensional spectra and by IFUs. Using Monte 
Carlo radiative transfer techniques, we study the nature 
of Mg II and Fe II absorption and emission for winds with 
a range of properties, accounting for the effects of reso- 
nant scattering and fluorescence. Although the winds are 
idealized, the results frequently contradict our intuition 
and challenge the straightforward conversion of observ- 
ables to (even crude) physical constraints. These findings 
have a direct bearing on recent and upcoming surveys of 
galactic outfiows, particularly those which make use of 
Mg II and Fe II transitions to probe outflow properties. 

The paper is organized as follows. In § [21 we describe 
the methodology of our radiative transfer algorithms. 
These are applied to a fiducial wind model in § [3] and 
variations of this model in § |4l In § [5l we explore wind 
models with a broader range of density and velocity laws. 
We discuss the principal results and connect to observa- 
tions in § [B] A brief summary is given in § [71 

2. METHODOLOGY 

This section describes our methodology for generating 
emission/absorption profiles from simple wind models. 

2.1. The Radiative Transitions 

In this paper, we focus on two sets of radiative transi- 
tions arising from Fe"*" and Mg+ ions (Table [TJ Figure [T]). 
This is a necessarily limited set, but the two ions and 
their transitions do have characteristics shared by the 
majority of low-ion transitions observed in cool-gas out- 
flows. Therefore, many of the results that follow may be 
generalized to observational studies that consider other 
atoms and ions tracing cool gas. 

The Mg"*" ion, with a single 3s electron in the ground- 
state, exhibits an alkali doublet of transitions at A ~ 
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Fig. 1.— Energy level diagrams for the Mg II AA2796, 2803 dou- 
blet and the UVl multip let of Fe II transitions (based on Figure 7 
from iHartigan et al.l 111999') V Each transition shown is labeled by 
its rest wavelength (A) and Einstein A-coefHcient (s~^). Black 
upward arrows indicate the resonance-line transitions, i.e. those 
connected to the ground state. The 2p®3p configuration of Mg+ is 
split into two energy levels that give rise to the Mg II AA2796, 2803 
doublet. Both the 3d® 4s ground state and 3d®4p upper level of Fc+ 
exhibit fine-structure splitting that gives rise to a series of electric- 
dipole transitions. The downward (green) arrows show the Fe II* 
transitions that arc connected to the resonance-line transitions (i.e. 
they share the same upper energy levels). We also show a pair of 
transitions (Fe II* AA2618, 2631; red and dashed lines) that arise 
from higher levels in the z ^D" configuration. These transitions 
have not yet been observed in galactic-scale outflows and are not 
considered in our analysis. 

2800A analogous to the Lya doublet of neutral hydro- 
gen. Figure [T] presents the energy level diagram for this 
Mg II AA2796, 2803 doublet. In non-relativistic quan- 
tum mechanics, the 2p^3p energy level is said to be split 
by spin-orbit coupling giving the observed line doublet. 
These are the only Mg II electric-dipole transitions with 
wavelengths near 2800A and the transition connecting 
the ^P3/2 and states is forbidden by several selec- 

tion rules. Therefore, an absorption from 2p^3s 2p^3p 
is followed w 100% of the time by a spontaneous decay 
(idecay ~ 4 X 10~^s) to the grouud state. Our treatment 
will ignore any other possibilities (e.g. absorption by a 
second photon when the electron is at the 2p^3p level). 

In terms of radiative transfer, the Mg II AA2796, 2803 
doublet is very similar to that for H I Lya, the 
Na I AA5891,5897 doublet, and many other doublets 
commonly studied in the interstellar medium (ISM) of 
distant galaxies. Each of these has the ground-state 
connected to a pair of electric dipole transitions with 
nearly identical energy. The doublets differ only in their 
rest wavelengths and the energy of the doublet separa- 
tion. For H I Lya, the separation is sufficiently small 
(Aw — cAE/E f« 1.3km s~^) that most radiative trans- 
fer treatments actually ignore it is a doublet. This is 
generally justifiable for Lya because most astrophysi- 
cal processes have turbulent motions that significantly 
exceed the doublet's velocity separation and effectively 
mix the two transitions. For Mg II (Aw « 770km s~^), 
Na I (Aw 304km s~^), and most of the other dou- 
blets commonly observed, the separation is large and the 
transitions must be treated separately. 



Iron exhibits the most complex set of energy levels for 
elements frequently studied in astrophysics. The_Fe+ 
ion a lone has millions of energy levels recorded (jKurucd 
l2005f l. and even this is an incomplete list. One reason 
for iron's complexity is that the majority of its config- 
urations exhibit fine-structure splitting. This includes 
the ground-state configuration (a ®D°) which is split 
into 5 levels, labeled by the total angular momentum 
J, with excitation energies Tex = AE/k ranging from 
Tex ~ 500-1500K (Figure[T]). Transitions between these 
fine-structure levels are forbidden (magnetic-dipole) and 
have spontaneous decay times of several hours. 

In this paper, we examine transitions between the 
ground-state configuration and the energy levels of the 
z configuration. This set of transitions (named the 
UVl multiplet) have wavelengths near 2600A. There are 
two resonance-line transitionf0 associated with this mul- 
tiplet (Fe II AA2586, 2600) corresponding to AJ = 0, -1; 
these are indicated by upward (black) arrows in Figure [TJ 
The solid (green) downward arrows in Figure [T] mark the 
non-resonant Fe II* transitions that are connected to the 
upper energy levels of the resonance lines. These transi- 
tions may occur following the absorption of a single pho- 
ton by Fe+ in its ground-state. This process may also be 
referred to as fluorescence. Note that two of these tran- 
sitions (Fe II* AA2626, 2632) are close enough in energy 
that their line profiles can overlap. 

The Figure also shows (as dashed, downward arrows) 
two of the Fe II* transitions that connect to higher en- 
ergy levels of the z configuration. Ignoring colli- 
sions and recombinations, these transitions may only oc- 
cur after the absorption of two photons: one to raise the 
electron from the ground-state to an excited state and 
another to raise the electron from the excited state to 
one of the z levels with J < 5/2. The excitation of 
fine-structure levels by the absorpt ion of UV photons is 
term ed indirect UV p umping (e.g ISilva fc Viegaal [200l 
Proc haska et al.l 120061 ) and requires the ion to lie near 
an intense source of UV photons. Even a bright, star- 
forming galaxy emits too few photons at A « 2500A to 
UV-pump Fe"*" ions that are farther than ~ 100 pc from 
the stars. In the following, we will assume that emission 
from this process is negligible. 

Our calculations also ignore collisional processetQ, i.e. 
coUisional excitation and de-excitation of the various 
levels. For the fine-structure levels of the a con- 
figuration, the excitation energies are modest (Tex ~ 
1000 K) but the critical density ric is large. For the 
a ®D°g/2 — ^a ^D°7/2 transition, the critical density w 
4 X lO^cm"'^. At these densities, one would predict de- 
tectable quantities of Fe I which has not yet been ob- 
served in galactic-scale outflows. If collisional excitation 
is insigniflcant then one may also neglect collisional de- 
excitation. Furthermore, observations rarely show ab- 
sorption from the flne-structure levels of the a config- 
uration and that material is not significantly blue-shifted 
(Rubin et al., in prep). In the following, we assume that 
electrons only occupy the ground-state, i.e. the gas has 
zero opacity to the non-resonant lines. Regarding the 

^ We adopt the standard convention that a "resonance line" is an 
electric-dipole transition connected to the ground-state. We also 
label non-resonant transitions with an asterisk, e.g. Fe II* A2612. 

® Recombination is also ignored. 
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Mg II doublet, its excitation energy is significantly higher 
implying negligible collisional processes at essentially any 
density. 

2.2. The Source 

Nearly all of the absorption studies of galactic-scale 
outflows have focused on intensely, star-forming galax- 
ies. The intrinsic emission of these galaxies is a com- 
plex combination of light from stars and H II regions 
that is then modulated by dust and gas within the ISM. 
For the spectral regions studied here, the hottest stars 
show a featureless continuum, but later spectral types 
do show significant Mg II and Fe II absorption. In addi- 
tion, asymmetric and/or blue-shifted absorption is exhib- 
ited in these transitions in A and F stars driving stellar 
winds (jSnow et al.lll994[ ). Mg II P Cygni profiles a re ob- 
served in a handful of F stars bv lSnow et al.l ()1994D . who 
attribute the emission to chromospheric activity rather 
than mass-loss effects. H II regions, meanwhile, are ob- 
served to emit weakly at the Mg II AA2796, 2803 dou- 
blet, primarily due to recombinations in the outer layers 
(|Kinnev et al.|[T993l ). It is beyond the scope of this paper 
to properly model the stellar absorption and H II region 
emission, but the reader should be aware that they can 
complicate the observed spectrum, independently of any 
outflow, especially at velocities v w km s~^. In the fol- 
lowing, we assume a simple flat continuum normalized 
to unit value. The size of the emitting region rgource is 
a free parameter, but we restrict its value to be smaller 
than the minimum radial extent of any gaseous compo- 
nent. Lastly, the source does not absorb any scattered 
or emitted photons. 

2.3. Monte Carlo Algorithm 

We calculated spectra using a 3D Monte Carlo 
radiation transport code originally d esigned for su- 
pernova outflows (jKasen et al.l |2006[) but modified 
to treat resonant line transport on a galactic scale 
(jKasen et al.l 120101 ). Our methods are similar to 
those used in several other codes which foc us on 
Lya (e.K.,lZ hcng & Miralda- Escudc 2002 ; Diikstr a et al.l 
l2006tlVerhamme et aLll2006HLaursen et al.ll2009l) . except 
that here we include the effects of multiple line scattering 
and fluorescence. 

For our 3-dimensional calculations, the wind properties 
(density, temperature and velocity) were discretized on 
a 300'^ Cartesian grid (for certain 1-D calculations, the 
wind properties were simply computed on the fly from 
analytic formulae). The radiation field was represented 
by N photon packets (typically N 10^) which were 
initially emitted isotropically and uniformly throughout 
the source region (r < Tsource)- The wavelengths of the 
packets were sampled from a flat spectral energy distri- 
bution. 

Photon packets were tracked through randomized scat- 
tering and absorption events until they escaped the com- 
putational domain. The distance to the next packet in- 
teraction event was determined by Doppler shifting pack- 
ets to the comoving frame and sampling the mean free 
path to each resonance line (and dust, if present). The 
resonant line opacity followed a Voig t profile wh i ch wa s 
determined using the analytic fits of iTasitsiomil ()2006[ ). 
Non-resonant lines were assumed to be completely opti- 
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Fig. 2. — Line profiles of a test problem calculated using the 
monte carlo code (solid) and by direct numerical integration of 
the transport equation under the Sobolev approximation (dashed 
lines). The configuration was spherical and consisted of a homolo- 
gously expanding wind with density profile n oc lying above an 
opaque emitting surface at velocity coordinate 100 km s~^. The 
wind extended to 2,000 km and the Sobolev optical depth had 
values of T = 1 (black lines), t = 10 (red lines), and r = 100 
(blue lines). Discrepancies in the r = 100 ca se are likely due to the 
breakdown of the Sobolev approximation (§ 15.1. Ill at high optical 
depth. 

cally thin, and the dust opacity was assumed to be wave- 
length independent and completely absorbing (§ 12. 4p . 
Because of the relatively low optical depths encountered 
in these models, no "core-skipping" scheme was applied 
to accelerate the resonant line transport, and recoil ef- 
fects were ignored. 

In a resonant line interaction, a photon excites an atom 
from the ground state to an excited level. The end re- 
sult is either scattering (i.e., de-excitation back to the 
ground state) or fluorescence (de-excitation to another 
excited level). The probability that the atom de-excites 
to lower level x is Pux = A^x/Y^iA^i, where the A^i 
are the Einstein spontaneous emission coefficients and 
the sum runs over all levels accessible from the upper 
state u. For each interaction event, a random number 
was drawn to sample the final state from this probability 
distribution. If the result was a scattering, the process 
was assumed to be coherent in the comoving frame. If 
the result was fiuorescence, the packet was reemitted at 
the line center wavelength of the new transition. Natu- 
ral line broadening in fiuorescence was ignored given the 
high velocity gradients in our models. In all cases, the 
angular redistribution was assumed to be isotropic. 

To generate multi-dimensional images and spectra of 
the system, an escape probability method was used. 
At the initiation of a packet, and at every subsequent 
interaction event, we calculated the probability Px — 
e^'^^'^Pux, that the packet de-excited into state x and es- 
caped the domain in some pre-specified direction. Here 
Tux is the optical depth to infinity, which was constructed 
by integration along the path to escape, taking into ac- 
count the relevant Doppler shifts and possible absorption 
by lines or dust. The contribution of the packet to the 
final spectrum was then added in for every possible final 
state, each shifted to the proper observer-frame wave- 
length and weighted by the probability Px- 
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To validate the monte carlo code, we calculated a series 
of test problems consisting of a point source embedded in 
a spherical, homogenous medium of varying resonant line 
optical depth. Th e resulting resonant line profiles (see 
IKasen et al.l[20T 0>) were found to be i n good agreement 
with the standard analytic solutio ns (jHarrington 1973,; 
iNeufeldl [T990I: iDiikstra et al.ll2006[ ). As a further test, 
we calculated line profiles for the case of an extended 
spherical source embedded in a homologously expanding 
wind with large velocity gradient. These model spectra 
(Figure [2|) are in good agreement with ones determined 
by direct numerical integration of the radiation transport 
equation under the Sob olev approximation fij 15.1.11 e.g., 
IJefferv fc BranchlfTOOOh . 

2.4. Dust 

For the majority of models studied in this paper, we 
assume the gas contains no dust. This is an invalid 
assumption, especially for material associated with the 
ISM of a galaxy. Essentially all astrophysical environ- 
ments that contain both cool gas and metals also show 
signatures of dust depletion and extinction. This in- 
cludes th e ISM of star- forming and H I-selected galax- 
ies ( e.g. ISavage fc SembachI [l99l iProchas ka fc Wolfd 
1200 it IProchaska eral...2007,l. strong Mg II metal-line ab - 



sorption systems (|York et al.ll2006riMenard etaIll2008D . 
and the galactic w i nds traced by low- ion transitions 
(|Pettini et all l2002t iR^pke et al.M2005b[ ). Extraplanar 
material likely associated with a galactic-scale outflow 
has been observed to emit I R radiation characteris- 
tic of dust (e. g. Hughej_eLalJ [1990; Ahon et al. 1999; 
IRadovich etaI.i,20Gl. ). In addition. .Menard et al., (,2G1Q1 
have argued from a statistical analysis that dust is dis- 
tributed to many tens of kpc from z ~ 0.1 galaxies and 
have suggested it was transported from the galaxies by 
galactic-scale winds. Although the galactic winds traced 
specifically by Mg II and Fe II transitions have not (yet) 
been demonstrated to contain dust, it is reasonable to 
consider its effects. 

For analysis on normalized spectra (i.e. absorption 
lines), the effects of dust are largely minimized; dust has 
a nearly constant opacity over small spectral regions and 
all features are simply scaled together. For scattered and 
resonantly trapped photons, however, the relative effect 
of dust extinction can be much greater. These photons 
travel a much longer distance to escape the medium and 
may experience a much higher integrated opacity from 
dust. Indeed, dust is frequently invoked to explain the 
wea k (or absent) Lya em ission from star-forming galaxies 
(e.g. lShaplev et all2003| ). Although the transitions stud- 
ied here have much lower opacity than Lya, dust could 
still play an important role in the predicted profiles. 

In a few models, we include absorption by dust under 
the following assumptions: (i) the dust opacity scales 
with the density of the gas (i.e. we adopt a fixed dust-to- 
gas ratio); (ii) the opacity is independent of wavelength, a 
reasonable approximation given the small spectral range 
analyzed; (iii) dust absorbs but does not scatter photons; 
(iv) the photons absorbed by dust are re-emitted at IR 
wavelengths and are 'lost' from the system. The dust 
absorption is normalized by Tdust , the integrated opacity 
of dust from the center of the system to infinity. The 
ambient ISM of a star-forming galaxy may be expected 
to exhibit Tdust values of one to a few at A '--^ 3000A (e.g. 



ICharlot fc FaUl[2000l) . 

2.5. The Sobolev Approximation 

A photon propagating through a differentially expand- 
ing medium interacts with a line only when its comoving 
frame wavelength is Doppler-shifted into near resonance 
with the line center rest wavelength. For a wind with 
a steep velocity law (i.e. a large gradient dv/dr) and/or 
a narrow intrinsic profile (i.e. a small Doppler parame- 
ter), the spatial extent of the region of resonance may be 
much smaller than the length scale of the wind itself. The 
interaction can then be co nsidered to occur at a point. 
In this case. lSobolevI (|1960[ ) introduced a formalism that 
gives the line optical depth r(r) at a given point in terms 
of the density and velocity gradient of the flow at that 
radius. For a wind in homologous expansion, this optical 
depth is independent of the direction of propagation and 
is given by 

Tsir)^KiX°\dv/dr\-' , (1) 



where 



Kf 



-hnt{r) 



(2) 



is the integrated line opacity, the line-center rest wave- 
length, ft the oscillator strength and the density in the 
lower level of the transition. We have neglected correc- 
tions for stimulated emission. This optical depth applies 
to a photon with wavelength A5 = ~ v{r) f ■ d/c) 
where d is the direction of propagation. The probability 
that such a photon is scattered/absorbed at the point 
of resonance is simply 1 — exp[— r5(r')]. We find that 
the Sobolev approximation applies for nearly all of the 
models presented in this paper, and therefore provides a 
convenient approach to estimating the optical depth. 

3. THE FIDUCIAL WIND MODEL 

In this section, we study a simple yet illustrative wind 
model for a galactic-scale outflow. The properties of this 
wind were tuned, in part, to yield a Mg II absorption 
profile si milar to those observed for z ^ 1 , star-fo rming 
galaxies (jWeiner et all 120091 : iRubin et al.l l2010a[ ). We 
emphasize, however, that we do not favor this fiducial 
model over any other wind scenario nor do its properties 
have special physical motivation. Its role is to establish 
a baseline for discussion. 

The fiducial wind is isotropic, dust-free, and extends 



from an inner wind radius ri^nc 
'"outer- It follows a density law. 



to an outer wind radius 



vor 



niiir) = r 

and a velocity law with a purely radial flow 

H — Vr{r)f 



''outer 



(3) 



(4) 



For the fiducial case, the hydrogen density at the inner 
radius is = 0.1 cm"^ and the velocity at the outer ra- 
dius is Wo — 1,000 km s~^. Turbulent motions are char- 
acterized by a Doppler parameteiQ foturb = 15km s~^. 



The adopted Doppler parameter has a minor impact on the 
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TABLE 2 
Wind Parameters: Fiducial Model 



Property 


Parameter 


Value 


TifinGTt"t/ 1 


n{r) 




Velocity law 


Vr 


oc r 


Inner Radius 


dinner 


1 kpc 


Outer Radius 


^^outcr 


20 kpc 


Source size 


^'source 


0.5 kpc 


Density Normalization 




0.1 Cm~'^ at Tinner 


Velocity Normalization 


vo 


1000km s^^at router 


Turbulence 


f'turb 


15 km s^-"- 


Mg+ Normalization 




10-5-*7nH 


Fe+ Normalization 




"Mg+/2 



We convert the hydrogen density njj to the number den- 
sities of Mg+ and Fe+ ions by assuming solar relative 
abundances with an absolute metallicity of 1/2 solar and 
depletion factors of 1/10 and 1/20 for Mg and Fe respec- 
tively, i.e. nMg+ — 10~^'^''nH and npg+=n-Mg+ /2. At the 
center of the wind is a homogeneous source of continuum 
photons with size ^source ■ The parameters for the fiducial 
wind model are summarized in Table [2j 

In Figure |3] we plot the density and velocity laws 
against radius; their simple power-law expressions are 
evident. The figure also shows the optical depth profile 
for the Mg II A2796 transition (T2796), estimated from the 
Sobolev approximatiorQ (Equation [1} . The T2796 profile 
peaks with r^gg ~ 30 at a velocity v,- w 65 km s~^ cor- 
responding to r w l.Srinner- The optical depth profile for 
the Mg II A2803 transition (not plotted) is scaled down 
by the /A ratio but is otherwise identical. Similarly, the 
optical depth profiles for the Fe II AA2586, 2600 transi- 
tions are scaled down by /A and the npg+/nMg+ ratio. 

For the spatially integrated (i.e., ID) spectra, the ac- 
tual dimensions and density of the wind are unimportant 
provided they scale together to give nearly the same op- 
tical depth. Therefore, one may consider the choices 
for nnner, ''outer, and as largely arbitrary. Never- 
theless, we adopted values for this fiducial model with 
some astrophysical motivation, e.g., values that corre- 
spond to galactic dimensions and a normalization that 
gives T^l"^ - 10. 

Using the methodology described in § [2l we propa- 
gated photons from the source and through the out- 
flow to an 'observer' at r ^ ''outer who views the entire 
wind-l-source complex. Figure H] presents the ID spec- 
trum that this observer would record, with the unatten- 
uated flux normalized to unit value. The Mg II doublet 
shows the canonical 'P-Cygni profile' that characterizes 
a continuum source embedded within an outflow. Strong 
absorption is evident at Sv = (A/Aq — 1) < —50 km in 
both transitions (with equivalent widths = 3.0A 

and ty^^"'^ — 1.3 A) and each shows emission at positive 
velocities. For an isotropic and dust-free model, the total 

results. The absorption profiles arc insensitive to its value and 
varying feturb only tends to modify the widths and modestly shift 
the centroids of emission lines. 

® We have verified this approximation holds by calculating T2796 
first as a function of velocity by summing the opacity for a series 
of discrete and small radial intervals between rinne^ and router- 
We then mapped T279e onto radius using the velocity law (Equa- 
tion |4]l. This optical depth profile is shown as the black dotted line 
in Figure |3] 




0.01 1 . I 

1 10 
Radius (kpc) 

Fig. 3. — Density (dashed; red), radial velocity (dotted; blue), 
and Mg II A2796 optical depth profiles (solid and dotted; black) 
for the fiducial wind model (see Table [2] for details). The density 
and velocity laws are simple and power-laws; these curves 
have been scaled for plotting convenience. The optical depth pro- 
file was calculated two ways: (i) using the Sobolev approximation 
(Equation [T] solid black curve) and (ii) summing the opacity at 
small and discrete radial intervals in velocity space and then con- 
verting to radius with the velocity law (doted black curve). These 
give very similar results. The wind parameters were set to give an 
optically thick medium at the inner radius (rinnor = 1 kpc) that 
becomes optically thin at the outer radius (router = 20 kpc). 



TABLE 3 
Equivalent Widths for the 
Fiducial Model 



Transition 


mt 

(km s^-"-) 


(4 


Mgll 2796 


[-1030, -65] 


2.96 




[-65,311] 


-1.75 


Mgll 2803 


[-458, -57] 


1.29 




[-57,692] 


-2.34 


Fell 2586 


[-365, -46] 


0.59 




[-46, 186] 


-0.10 


Fell 2600 


[-626, -78] 


1.16 




[-49, 470] 


-0.94 


Fell* 2612 


[-219,183] 


-0.33 


Fell* 2626 


[-194,177] 


-0.27 


Fell* 2632 


[-155, 130] 


-0.16 



^ Velocity interval over which the 
equivalent width is calculated. 
^ Negative values indicate line- 
emission. 



equivalent width of the doublet must be zero, i.e. every 
photon absorbed eventually escapes the system, typically 
at lower energy. The wind simply shuffles the photons in 
frequency space. A simple summation of the absorption 
and emission equivalent widths (Table confirms this 
expectation. 

Focusing further on the Mg II absorption, one notes 
that the profiles lie well above zero intensity and have 
similar depth even though their /A values differ by a fac- 
tor of two. In standard absorption-line analysis, this is 
the tell-tale signature of a 'cloud' that has a high optical 
depth (i.e. is saturat ed) which on ly partially cov ers the 
emitting source (e.g. iRupke et al.,, 2005a; .Haniann et al.l 
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Wavelength (Ang) Relative Velocity (km s"') 



Fig. 4.— Left - (Upper) Mg II AA2796, 2803 profiles for the fiducial wind model described in Tableland Figure |3] The doublet shows 
the P-Cygni profiles characteristic of an outflow with significant absorption blueward of line-center (dashed vertical lines) extending to 
V = —1000km s~^ and significant emission redward of line-center. Note that even though the peak optical depth of the Mg II transitions 
is T^gg 30 at D « —70km s~^, photons scattered off the outflow significantly fill-in the absorption. (Lower) Fe II absorption and 
emission profiles for the UVl multiplet at A R:i 2600A. The Fe II AA2586, 2600 resonance lines show weaker absorption due to the smaller 
Fe+ number density and lower /A values. Each also shows a P-Cygni profile, although the emission for Fe II A2586 is much weaker than 
that for the Fe II A2600 and Mg II AA2796, 2803 transitions. This is because a majority of the absorbed Fe II A2586 photons fluoresce into 
Fe II* AA2612, 2632 photons. Right - A subset of the transitions displayed in a velocity plot. 



[2OT0l ). Our fiducial wind model, however, entirely cov- 
ers the source; the apparent partial covering must be 
related to a different effect. Figure [5] further empha- 
sizes this point by comparing the absorption profiles from 
Figure |4] against an artificial model where no absorbed 
photons are re-emitted. As expected from the T2796 pro- 
file (Figure 13]), this 'intrinsic' model produces a strong 
Mg II AA2796,2803 doublet that absorbs all photons 
at Sv w -100 km s~\ i.e. I^^q = exp(-T2"^|g) w 0. 
The true model, in contrast, has been 'filled in' at 
Sv w —100 km s""'^ by photons scattered in the wind. An 
absorption-line analysis that ignores these effects would 
(i) systematically underestimate the true optical depth 
and/or (ii) falsely conclude that the wind partially cov- 
ers the source. We will find that these are generic results, 
even for wind models that include dust and are not fully 
isotropic. 

Turning to the emission profiles of the 
Mg II AA2796,2803 doublet, one notes that they 
are quite similar with comparable equivalent widths. 



This is because the gas is optically thick yielding com- 
parable total absorption. The flux of the Mg II A2803 
transition even exceeds that for Mg II A2796 because the 
wind speed is greater than the velocity separation of the 
doublet, |wr|max > (A'y)Mgii- Therefore, the red wing of 
the Mg II A2796 emission profile is partially absorbed 
by Mg II A2803 and re-emitted at lower frequency. This 
yields a line ratio that is far below the 2 : 1 ratio that 
one may have naively expected (e.g. if the line-emission 
resulted from recombinations), and leads us to conclude 
that the relative strengths of the emission lines are 
sensitive to both the opacity and velocity extent of the 
wind. 

Now consider the Fe II transitions: the bottom left 
panel of Figure |3] covers the majority of the Fe II UVl 
transitions and several are shown in the velocity plot. 
The line profile for Fe II A2600 is very similar to the 
Mg II AA2796, 2803 doublet; one observes strong ab- 
sorption to negative velocities and strong emission at 
Sv > km s^^ producing a characteristic P-Cygni pro- 
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Fig. 5. — The solid curves show the line profiles (absorption and 
emission) of the Mg 11 and Fe II resonance lines for the fiducial v^ind 
model. These include the effects of scattered photons and show the 
canonical P-Cygni profiles of a source embedded within an outflow. 
Overplotted on each transition (dotted red line) is the predicted 
absorption profile under the constraint that every absorbed photon 
is lost from the system (referred to as the intrinsic profile), i.e. no 
scattering or re-emission occurs. The intrinsic profiles have been 
'fiUed-in' significantly with photons scattered by the wind into our 
sightline. Ignoring this process, one would model the absorption 
lines with a systematically lower optical depth and/or (incorrectly) 
conclude that the source is partly covered by the gas. 

file. Splitting the profile at dv = —50 km s^^, we measure 
an equivalent width W^^^'^ = 1.16 A in absorption and 
^^2600 ^ _o.94A in emission (Tabled for a total equiva- 
lent width of M^TOT ~ O.22A. In contrast, the Fe II A2586 
resonance line shows much weaker emission and a much 
higher total equivalent width (W|^^| = 0.49 A), even 
though the line has a 2x lower /A value. These dif- 
ferences between the Fe II resonance lines (and between 
Fe II and Mg II) occur because of the complex of non- 
resonant Fe II* transitions that are coupled to the res- 
onance lines (Figured]). Specifically, a resonance pho- 
ton absorbed at Fe II AA2586, 2600 has a finite prob- 
ability of being re-emitted as a non-resonant photon 
which then escapes the system without further interac- 
tion. The principal effects are to reduce the line emission 
of Fe II AA2586, 2600 and to produce non-resonant line- 
emission (e.g. Fe II* A2612). 

The reduced Fe II A2586 emission relative to 
Fe II A2600 is related to two factors: (i) there is an addi- 
tional downward transition from the z level and 
(ii) the Einstein A coefficients of the non-resonant lines 
coupled to Fe II A2586 are comparable to and even ex- 
ceed the Einstein A coefficient of the resonant transition. 
In contrast, the Fe II* A2626 transition (associated with 




-200 

V (km s"') 




10 20 -20 -10 

kpc kpc 
Fig. 6. — Surface-brightness emission maps around the 
Mg 11 A2796 transition for the source+wind complex of the fidu- 
cial model. The middle panel shows the ID spectrum with v = 
0km s~^ corresponding to A = 2796. 35A and the dotted verti- 
cal curves indicate the velocity slices for the emission maps. The 
source has a size rsource = 0.5 kpc, traced by a few pixels at the 
center of each map. At i; = —250 km , the wind has an optical 
depth of T2796 ~ 1 and the source contributes roughly half of the 
observed fiux. At v = —100km s~^ the wind absorbs all photons 
from the source and the observed emission is entirely due to pho- 
tons scattered by the wind. Amazingly, this emission exceeds the 
integrated flux at ii = -250km s~^ such that the Mg 11 A2796 
line center is offset from the (intrinsic) peak in optical depth. At 
V > km s~^, both the source and wind contribute to the observed 
emission. At all velocities, the majority of emission comes from the 
inner 5 kpc. 



Fe II A2600) has an approximately 4x smaller A coeffi- 
cient than the corresponding resonance line. Therefore, 
the majority of photons absorbed at A « 2600A are re- 
emitted as Fe II A2600 photons, whereas the majority of 
photons absorbed at A « 2586A are re-emitted at longer 
wavelengths (Fe II* A2612 or A2632). If we increase T2600 
(and especially if we include gas with w km s~^) then 
the Fe II A2600 emission is significantly suppressed (e.g. 
ij 14.31) . The total equivalent width, however, of the three 



lines connected to the z upper level must still 

vanish (photons are conserved in this isotropic, dust-free 
model). 

The preceding discussion emphasizes the filling-in of 
resonance absorption at 6v < —50 km s~^ and the gen- 
eration of emission lines at Sv « km s""'^ by photons 
scattered in the wind. We also mapped the emission of 
the fiducial model to sudy its spatial extent (see ii l2.3l for 
a description of the algorithm). The output is a set of 
surface-brightness maps in a series of frequency channels 
yielding a dataset analogous to integral-field-unit (IFU) 
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Fig. 7. — (Upper) Surface-brightness emission maps around the 
Fe II A2600 transition for the source+wind complex. The results 
are very similar to those observed for the Mg II AA2796, 2803 dou- 
blet (Figure I6J. (Lower) Surface-brightness emission maps around 
the Fe II A2612 transition for the source+wind complex. In this 
case, the source is unattenuated yet scattered photons from the 
wind also make a significant contribution. 

observations. In Figure [51 we present the output at sev- 
eral velocities relative to the Mg II A2796 transition. At 
Sv = —250km s""'^, where the wind has an optical depth 
T2796 < 1 (Figure [S]), the source contributes roughly half 
of the observed flux. At Sv = — 100km s~^, however, 
the wind absorbs all photons from the source and the 
observed emission is entirely from photons scattered by 
the wind. This scattered emission actually exceeds the 
source-|-wind emission at Sv — — 250kms~^ such that 
the absorption profile is negatively offset from the veloc- 
ity where T2796 is maximal (Figure [3]) . The net result 
is weaker Mg II absorption that peaks blueward of the 
actual peak in the optical depth profile. Table |4] reports 
several kinematic measurements of the absorption and 
emission features. Clearly, these effects complicate esti- 
mates for the speed, covering fraction, and total column 
density of the wind. At Sv = km s~^, the wind and 
source have comparable total flux with the latter domi- 
nating at higher velocities. 

Similar results are observed for the Fe II resonance 
transitions (Figure [7]) . For transitions to fine-structure 
levels of the a configuration, the source is unattenu- 
ated but there is a significant contribution from photons 
generated in the wind. 

At all velocities, the majority of light comes from the 
inner regions of the wind. The majority of Mg II emission 
occurs within the inner few kpc, e.g. 50 — 60% of the light 
at Sv = —100 to +100 km s^^ comes from |r| < 3kpc. 
The emission is even more centrally concentrated for the 
Fe II transitions. A proper treatment of these distri- 
butions is critical for interpreting observations acquired 
through a slit, i.e. where the aperture has a limited ex- 
tent in one or more dimensions. A standard longslit on 
lOm-class telescopes, for example, subtends w 1" corre- 
sponding to 5 — 10 kpc for z ~ 1. We return to this issue 
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Fig. 8. — Profiles of the Fe II and Mg II profiles for the fiducial 
case (black lines) compared against an anisotropic wind blowing 
into only 2it steradians as viewed from ij> = 0° (source uncovered) 
to (/i = 180° (source covered). One detects significant emission for 
all orientations but significant absorption only for (p > 90°. The 
velocity centroid of the emission shifts from positive to negative 
velocities as increases and one transitions from viewing the wind 
as lying behind the source to in front of it. The velocity centroid 
of emission, therefore, diagnoses the degree of anisotropy for the 
wind. 

in §[6l 

The results presented in Figures[n]and[7]are sensitive to 
the radial extent, morphology, density and velocity pro- 
files of this galactic-scale wind. Consequently, IFU ob- 
servations of line emission from low-ion transitions may 
offer the most direct constraints on galactic-scale wind 
properties. 

4. VARIATIONS TO THE FIDUCIAL MODEL 

In this section, we investigate a series of more com- 
plex wind scenarios through modifications to the fidu- 
cial model. These include relaxing the assumption of 
isotropy, introducing dust, adding an ISM component 
within rinnorj and varying the normalization of the opti- 
cal depth profiles. 

4.1. Anisotropic Winds 

The fiducial model assumes an isotropic wind with 
only radial variations in velocity and density. Angu- 
lar isotropy is obviously an idealized case, but it is fre- 
quently assumed in studies of galactic-scale outflows (e.g. 
iSteidel et al.ll2010[ ). There are several reasons, however, 
to consider anisotropic winds. Firstly, galaxies are not 
spherically symmetric; the sources driving the wind (e.g. 
supernovae, AGN) are very unlikely to be spherically dis- 
tributed within the galaxy. Secondly, the galactic ISM 
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Fig. 9. — Profiles of the Fe II and Mg II profiles for the fiducial 
case (black lines) compared against a bi-conical wind with varying 
opening angle {9i, = 0° is a full wind). For modest angles, one 
observes similar results to the fiducial model but larger angles {Bi, > 
45°) the emission is significantly suppressed at all velocities. 

frequently has a disk-like morphology which will suppress 
the wind preferentially at low galactic latitudes , per- 
haps yielding a bi-conic m orphology fe-g. ^Heckman et alJ 
119901 iWalter et aLll2002| ). Lastly, the galaxy could be 
surrounded by an aspherical gaseous halo whose interac- 
tion would produce an irregular outflow. 

With these considerations in mind, we reanalyzed the 
fiducial model with the 3D algorithm after departing 
from isotropy. It is beyond the scope of this paper to 
explore a full suite of anisotropic profiles. We consider 
two simple examples: (i) half the fiducial model, where 
the wind density is set to zero for 27r steradians (i.e. a 
hemispherical wind). This model is viewed from (/) — 0° 
(source uncovered) to = 180° (source covered); and (ii) 
a bi-conical wind which fills \9\ < Ojj both into and out of 
the plane of the sky and that is viewed along the axis of 
rotational symmetry {6 = 0° is defined along this axis). 

The resulting Mg II and Fe II profiles for the half wind 
are compared against the fiducial model (isotropic wind) 
in Figure m Examining the Mg II AA2796, 2803 doublet, 
the (/) = 0° model only shows line-emission from photons 
scattered off the back side. These photons, by definition, 
have 5v > km s^^ relative to line-center (a subset have 
Sv < km s~^ because of turbulent motions in the wind) . 
When viewed from the opposite direction (0 — 180°), the 
absorption lines dominate, but there is still significant 
line-emission at dv ~ km s~^ and at Sv < km s~^ 
from photons that scatter through the wind which fills 
in the absorption. The key difference between this and 
the isotropic wind is the absence of photons scattered to 



Sv > 100km s~^; this also implies deeper Mg II A2803 
absorption at 6v « —100km s~^. The shifts in velocity 
centroid and asymmetry of the emission lines serve to 
diagnose the degree of wind isotropy, especially in con- 
junction with analysis of the absorption profiles. The 
results are similar for the Fe II AA2586, 2600 resonance 
lines. The Fe II* A2612 line, meanwhile, shows most 
clearly the offset in velocity between the source unob- 
scured (</> — 0°) and source covered (0 — 180°) cases. 
The offset of the Fe II* lines is the most significant dif- 
ference from the fully isotropic wind. 

We have also analyzed the predicted profiles for a series 
of bi-conical winds with 9b = 30 — 80°, each viewed along 
the axis of rotational asymmetry with the wind covering 
the source. Figure |9] and Table |4] summarizes the re- 
sults for several cases. For 9i, = 80°, the results are very 
similar to the fiducial model; the profiles show similar 
equivalent widths for absorption and emission. The key 
quantitative difference is that the Fe II* line-emission is 
modestly suppressed at Sv w km s~^ because the por- 
tion of the wind with that projected velocity has been 
removed. For more highly collimated winds 9b < 45°, 
however, the line-emission is several times weaker than 
the fiducial model. Similarly, the absorption-line profiles 
more closely track the intrinsic optical depth of the out- 
flow. Of course, these same models when viewed 'side- 
ways' would show no absorption but strong line-emission. 
Such events are sufficiently rare (Rubin et al., in prep) 
that we consider this level of anisotropy to be uncom- 
mon. Nevertheless, the results for a bi-conic clearly have 
important implications for the nature of absorption and 
line-emission in our fiducial model. 

4.2. Dust 

As described in § 12.41 one generally expects dust in as- 
trophysical environments that contain cool gas and met- 
als. This dust modifies the observed wind profiles in 
two ways. First, it is a source of opacity for all of the 
photons. This suppresses the flux at all wavelengths by 
w exp(— Tdust) but because we re-normalize the profiles, 
this effect is essentially ignored. Second, photons that 
are scattered by the wind must travel a greater distance 
to escape and therefore suffer from greater extinction. 
A photon that is trapped for many scatterings has an 
increased probability of being absorbed by dust. Sec- 
tion 12.41 describes the details of our treatment of dust; 
we remind the reader here that we assume a constant 
dust-to-gas ratio that is normalized by the total opti- 
cal depth Tdust photons would experience if they traveled 
from the source to infinity without scattering. 

In Figure [TOl we show the Mg II and Fe II profiles of the 
fiducial model (rdust = 0) against a series of models with 
Tdust > 0. For the Mg II transitions, the dominant ef- 
fect is the suppression of line emission at Sv > km s~^. 
These 'red' photons have scattered off the backside of the 
wind and must travel a longer path than other photons. 
Dust leads to a differential reddening that increases with 
velocity relative to line-center. This is a natural con- 
sequence of dust extinction and is most evident in the 
Fe II* A2612 emission profile which is symmetrically dis- 
tributed around Sv = 0km in the Tdust = model. 
The degree of suppression of the line-emission is rela- 
tively modest, however. Specifically, we find that the 
flux is reduced by a factor of the order of (1+Tdust)~^ 
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Fig. 10. — Profiles of the Fe II and Mg II profiles for the fiducial 
model (black) compared against a series of models that include 
dust extinction, parameterize by the integrated dust opacity from 
source to infinity (rdugt). The primary effect of dust is to suppress 
the line-emission relative to the continuum. A more subtle but 
important effect is that the redder photons in the emission lines 
(corresponding to positive velocity offsets relative to line-center) 
suffer greater extinction. This occurs because the 'redder' photons 
that we view have travelled farther to scatter off the backside of the 
wind. Note that the absorption lines are nearly unmodified until 
"""dust = 10, a level of extinction that would preclude observing the 
source altogether. 

(Figure fTT|) . instead of the factor exp(-Tdust) that one 
may have naively predicted. In terms of absorption, the 
profiles are nearly identical for rjust < 3. One requires 
very high extinction to produce a deepening of the pro- 
files at Sv « —100 km s~^. 

We conclude that dust has only a modest infiuence on 
this fiducial model and, by inference, on models with 
moderate peak optical depths and significant velocity 
gradients with radius (i.e. scenarios in which the pho- 
tons scatter only one to a few times before exiting). For 
qualitative changes, one requires an extreme level of ex- 
tinction (rdust = 10). In this case, the source would 
be extinguished by 15 magnitudes and could never be 
observed. Even Tdust = 3 is larger than typically in- 
ferre d for the star-formin g galaxies that drive outflows 
(e.g. ICharlot fc FaUl[2000l) . For the emission lines, the 
dominant effect is a reduction in the flux with a greater 
extinction at higher velocities relative to line-center. In 
these respects, dust extinction crudely mimics the be- 
havior of the anisotropic wind described in Section 14.11 

4.3. ISM 

The fiducial model does not include gas associated 
with the ISM of the galaxy, i.e. material at r Okpc 
with Uj. w km s~^. This allowed us to focus on re- 



< 




Fig. 11.— Predicted equivalent widths for Mg II A2803 line- 
emission from the fiducial wind model, as attenuated by dust with 
a range of Tjugt values. The black (solid) curve shows the model 
results. The dashed (green) curve shows Wc{Td) = Wc{Td = 0) ■ 
exp(— Tjj). ft is obvious that the equivalent widths (i.e., the fiux 
relative to the continuum) do not follow this scaling. The dotted 
(red) line, meanwhile, plots We(Trf) = We{Td = 0) • (1 -I- r^j)"^. 
This simple approximation is a good representation of the results 
for our radiative transfer calculations. 



suits related solely to a wind component. The decision 
to ignore the ISM was also motivated by the general 
absence of significant absorption at Sv ^ km s~^ in 
galaxies that exhib i t outflows (e . g. iWeiner et al.l 120091 : 
IRubin et alJ I2010bl: ISteidel et al. l l2010D . On the other 
hand, the stars that comprise the source are very likely 
embedded within and fueled by gas from the ISM. Con- 
sider, then, a modiflcation to the fiducial model that in- 
cludes an ISM. Specifically, we assume the ISM compo- 
nent has density tih = 1 cm~^ for risM < r < r-mnci- with 
nsM — 0.5 kpc, an average velocity of v^- = 0km s^^, 
and a larger turbulent velocity 6ism = 40 km s^^. The 
resultant optical depth profile T2796 is identical to the 
fiducial model for r > 2 kpc, has a slightly higher opac- 
ity at r = 1 — 2 kpc, and has a very large opacity at 
r = 0.5 — 1 kpc. 

In Figure [121 the solid curves show the Mg II and Fe II 
profiles for the ISM-|-wind and fiducial models. In com- 
parison, the dotted curve shows the intrinsic absorption 
profile for the ISM-t-wind model corresponding to the 
(unphysical) case where none of the absorbed photons are 
scattered or re-emitted. Focus first on the Mg II doublet. 
As expected, the dotted curve shows strong absorption 
at Sv ~ km and blueward. The full models, in con- 
trast, show non-zero fiux at these velocities and even a 
normalized flux exceeding unity at Sv w km s~^. In 
fact, the ISM-t-wind model is nearly identical to the fidu- 
cial model; the only quantitative difference is that the 
velocity centroids of the emission lines are shifted red- 
ward by « -hlOOkms^^. We have also examined the 
spatial distribution of emission from this model and find 
results qualitatively similar to the fiducial wind model 
(Figures Eland [ZD. 

There are, however, several qualitative differences for 
the Fe II transitions. First, the Fe II A2586 transition 
in the ISM-fwind model shows much stronger absorp- 
tion at (5u ~ to —100 km s~^. In contrast to the Mg II 
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Fig. 12.— Profiles of the Fe II and Mg II profiles for the 
ISM+wind model (red) compared against the fiducial wfind model 
(black). The dotted line traces the predicted absorption pro- 
file in the absence of re-emission and scattering. Regarding the 
Mg II doublet, the primary difference between the ISM-|-wind and 
the fiducial models is the shift of f» 100 km in the emission lines 
from V Si km s~^ for the fiducial model to v ^ -1-100 km s~^ for 
the ISM-fwind model. The Fe II profiles, meanwhile, show several 
qualitative differences. The Fe II AA2586, 2600 resonance transi- 
tions each exhibit much greater absorption a,t v km than 
the fiducial model. The resonant line-emission is also substantially 
reduced, implying much higher fiuxes for the non-resonant lines 
(e.g. Fe II* A2612,2626). Lastly, we note that the Fe II A2586 ab- 
sorption profile provides a good (albeit imperfect) representation 
of the wind opacity and therefore offers the best characterization 
of an ISM component. 

doublet, the profile is not filled in by scattered photons. 
Instead, the majority of Fe II A2586 photons that are ab- 
sorbed are re-emitted as Fe II*AA2612, 2632 photons. In 
fact, the Fe II A2586 profile very nearly matches the pro- 
file without re-emission (compare to the dotted lines); 
this transition provides a very good description of the 
intrinsic ISM-|-wind optical depth profile. We conclude 
that resonant transitions that are coupled to (multiple) 
non-resonant, electric dipole transitions offer the best di- 
agnosis of ISM absorption. 

The differences in the Fe II absorption profiles are re- 
flected in the much higher strengths (3 — lOx) of emis- 
sion from transitions to the excited states of the a 
configuration. This occurs because: (1) there is greater 
absorption by the Fe II AA2586, 2600 resonance lines; 
and (2) the high opacity of the ISM component leads 
to an enhanced conversion of resonance photons with 
Sv « km into Fe II* photons. This is especially no- 
table for the Fe II A2600 transition whose coupled Fe II* 
transition shows an equivalent width nearly 10 x stronger 
than for the fiducial model. The relative strengths of the 



Fe II A2600 and Fe II* A2626 lines provide a direct diag- 
nostic of the degree to which the resonance line photons 
are trapped, i.e. the peak optical depth of Fe II 2600 and 
the velocity gradient of the wind. 

Because of the high degree of photon trapping within 
the ISM component, this model does suffer more from 
dust extinction than the fiducial model. We have stud- 
ied the ISM-|-wind model including dust with Tdust = 1 
(normalized to include the ISM gas). All of the emission 
lines are significantly reduced. The Mg II emission is af- 
fected most because these lines are resonantly trapped. 
The Fe II* emission is also reduced relative to the dust- 
free model, but the absolute fiux still exceeds the fiducial 
model (Table [4]). At large negative velocity offsets from 
systemic, the two profiles are nearly identical. We con- 
clude that a dusty ISM model could show significant ab- 
sorption at Sv ~ km s~^ in resonance lines (including 
Mg II) with strong line-emission in the Fe II* transitions. 
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Fig. 13. — Mg II and Fe II profiles for the fiducial model with 
varying normalization, parameterized by n^. As expected, the 
strength of absorption increases with increasing n%; this also re- 
sults in stronger line-emission. Note that the Fe II A2586 emission 
is always weak. Only its absorption increases with and actually 

exceeds the depth of Fe II A2600 for n° > n°'^. The depth of the 
Mg II doublet, meanwhile, always falls below a relative fiux of 0.3 
while the Mg II A2803 emission rises steadily with n^. 



4.4. Varying 

The final modification to the fiducial model considered 
was a uniform variation of the normalization of the op- 
tical depth profiles. Specifically, we ran a series of addi- 
tional models with at 1/3, 3, and 10 times the fiducial 



value of n 



oj 



0.1 cm ^. The resulting Mg II and Fe II 
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Line Diagnostics for the Fiducial Model and Variants 
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Note. — Listed are the equivalent widths (intrinsic, absorption, and emission), the peak optical depth for the absorption Tpk = — ln(/niin), the velocity where 
the optical depth peaks v-j-, the optical depth-weighted velocity centroid = f dvv\n[I{v)]/ f dv\n[I(v)], the peak flux /pk in emission, the velocity where the 
flux peaks v/, the flux- weighted velocity centroid of the emission line (occasionally affected by blends with neighboring emission lines), and the 90% width 
Ave- The and vf^^ columns give the velocity range used to calculate the absorption and emission characteristics, respecitvely. These were defined by the 
velcoities where the profile crossed 0.95 in the normalized fiux. 
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profiles are compared against the fiducial model in Fig- 
ure [TS] Inspecting Mg II, one notes that the models with 
1/3 and 3n^-^ behave as expected. Higher/lower optical 
depths lead to greater/ weaker Mg II A2796 absorption 
and stronger/weaker Mg II A2803 emission. In the ex- 
treme case of = lOn^^ , the Mg II A2803 absorption 
and Mg II A2796 emission have nearly disappeared and 
one primarily observes very strong Mg II A2796 absorp- 
tion and Mg II A2803 emission. In essence, this wind 
has converted all of the photons absorbed by the Mg II 
doublet into Mg II A2803 emission. 

Similar behavior is observed for the Fe II transitions. 
Interestingly, the Fe II A2586 transitions never show sig- 
nificant emission, only stronger absorption with increas- 
ing rij\. In fact, for nj\ > 2^^''^ the peak optical depth of 
Fe II A2586 actually exceeds that for Fe II A2600 because 
the latter remains significantly filled-in by scattered pho- 
tons. By the same token, the fiux of the Fe II* A2612 
emission increases with n^. 

4.5. Summary Table 

Table |4] presents a series of quantitative measures of 
the Mg II and Fe II absorption and emission lines for the 
fiducial model (§[3]) and a subset of the models presented 
in this section. Listed are the absorption and emission 
equivalent widths (Wa,We), the peak optical depth for 
the absorption Tpk, the velocity where the optical depth 
peaks Vr, the optical depth- weighted velocity centroid 
Vf = J dvvhi[I{v)]/ J dv\Ti[I{v)], the peak fiux /pk in 
emission, the velocity where the fiux peaks w/, and the 
fiux- weighted velocity centroid of the emission line vj. 
We discuss several of these measures in § El 

5. ALTERNATE WIND MODELS 

This section presents several additional wind scenarios 
that differ significantly from the fiducial model explored 
in the previous sections. In each case, we maintain the 
simple assumption of isotropy. 

5.1. The Lyman Break Galaxy Model 

The Lyman break galaxies (LBGs), UV color-selected 
galaxies at z ~ 3, exhibit cool gas outflows in Si II, 
C II, etc. transitions with speeds up to 1000 km s~^ 
(e.g. iLowenthal et all 119971: iPettini et all I1998D . Re- 
searchers have invoked these winds to explain enrichment 
of the intergalactic me dium (e.g. .A guirre e t al. 2001]; 
IScannapieco et all 120061), the origin of the damped Lya 
svstems ( |Nulsen et al.l ll998HSchavel l200lD. and the for- 
matio n of 'red and dead' galaxies Ce.g. iSomerville et al.l 
I2001D . Although the presence of these outflows was es- 
tablished over a decade ago, the processes that drive 
them remain unidentifled. Similarly, current estimates 
of the mass and energetics of the outflows suffer from 
orders of m agnitude uncer t ainty. 

Recently, iSteidel et afl (|2010l . hereafter SIO) intro- 
duced a model to explain jointly the average absorp- 
tion they observed along the sightlines to several hun- 
dred LBGs and the average absorption in gas observed 
transverse to these galaxies. Their wind model is deflned 
by two expressions: (i) a radial velocity law v,-{r); and 
(ii) the covering fraction of optically thick 'clouds' fc{r)- 
For the latter, SIO envision an ensemble of small, opti- 
cally thick (r » 1) clouds that only partially cover the 



galaxies. For the velocity law, they adopted the following 
functional form: 



^LBG 

I — a 



1/2 



,1/2 



(5) 



with Albg the constant that sets the terminal speed, 
Hnncr the inner radius of the wind (taken to be 1 kpc) , 
and a the power-law exponent that describes how steeply 
the velocity curve rises. Their analysis of the LBG ab- 
sorption proflles implied a very steeply rising curve with 
a w 1.3. This velocity expression is shown as a dotted 
line in Figure [Tik . 

The covering fraction of optically thick clouds, mean- 
while, was assumed to have the functional form 



(6) 



with 7 0.5 and fc,max the maximum covering frac- 
tion. From this expression and the velocity law, one can 
recover an absorption profile /lbg(''^) — l~/c[?'(i')], writ- 
ten explicitly as 



'LBG 



iv)^l- fc. 



1 



A 



LBG 



7/(a-l) 



(7) 

The resulting profile for fc,max = 0.6, 7 = 0.5, a — 
1.3, and Albg = -192, 000 km^ s^^ kpc^^ is displayed 
in Figure [Hb. 

In the following, we consider two methods for analyz- 
ing the LBG wind. Both approaches assume an isotropic 
wind and adopt the velocity law given by Equation [5l In 
the first model, we treat the cool g diffuse medium 

with unit covering fraction and a radial density profile 
determined from the Sobolev approximation. We then 
apply the Monte Carlo methodology used for the other 
wind models to predict Mg II and Fe II line profiles. 
In the second LBG model, we modify our algorithms to 
more precisely mimic the concept of an ensemble of op- 
tically thick clouds with a partial covering fraction on 
galactic scales. 

5.1.1. Sobolev Inversion 

As demonstrated in Figure [141 the wind velocity for 
the LBG model rises very steeply with increasing radius 
before fiattening at large radii. Under these conditions, 
the Sobolev approximation provides an accurate descrip- 
tion of the optical depth for the fiowing medium. The 
Sobolev line optical depth profile can be derived from the 
absorption profile (Equation [7]), 



TLBG (t') = - In [ -^LBG {v) ] 



(8) 



This simple expression assumes that the source size is 
small compared the wind dimensions and neglects the 
effects of light scattered in the wind (for more general in- 
versi on formulae th at do not make these assumptions, see 
Kase n et al.l (|2002l )). The Sobolev approximation (Equa- 
tion [T|) provides a simple expression for the optical depth 
at each radius in terms of the local density and velocity 
gradient. In general, the velocity gradient must be taken 
along the direction of propagation of a photon, however 
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Fig. 14. — (a) The dotted curve shows the velocity law (la- 
beled on the right axis) for the LBG model of SIO plotted against 
(r/kpc — 1) to highlight the evolution in the quantities at small 
radii. For example, note how rapidly the velocity rises from r = 1 
to 2 kpc. The solid curve shows the Sobolev solution for the Mg+ 
gas density derived from the average absorption profile of LBG 
galaxies (see below). The density falls off initially as (r/kpc— l)""'^ 
and then steepens to (r/kpc — l)"'^'*. (b) The average profile for 
LBG cool gas absorption as measured and defined by SIO (black 
solid curve). Overplotted on this curve is a (red) dotted line that 
shows the Mg II A2796 absorption profile derived from the density 
(and velocity) law shown in the upper panel. 

in this model the central source is relatively small and 
the radiation field is therefore primarily radially directed. 
In this case the velocity gradient in the radial direction, 
dv-c/ dr, is the relevant quantity, and we may use Equa- 
tion [I] to invert the optical depth and determine the den- 
sity profile. 



\dv^/dr\ I rrieC 
ULBGir) = tlbg(0 j^^o [ ^ 



(9) 



The solid curve in Figure [Tik shows the resultant den- 
sity profile for Mg+ assuming that the Mg II A2796 line 
follows the intensity profile drawn in Figure [Mb . This is 
a relatively extreme density profile. From the inner ra- 
dius of 1 kpc to 2 kpc, the density drops by over 2 orders 
of magnitude including nearly one order of magnitude 
over the first 10 pc. Beyond 2 kpc, the density drops 
even more rapidly, falling orders of magnitude from 2 to 
100 kpc. 

We verified that the Sobolev-derived density profile 
shown in Figure 114b reproduces the proper absorption 
profile by discretizing the wind into a series of layers and 
calculating the integrated absorption profile. This calcu- 
lation is shown as a dotted red curve in Figure [Tib : it is 
an excellent match to the desired profile (black curve). 
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Fig. 15.— Profiles of the Fe II and Mg II profiles for the two LBG 
models considered: (red) the profiles predicted using the Sobolev 
approximation and (black) a model constructed to faithfully rep- 
resent the wind model advocated by SIO. The dotted line shows 
the absorption profiles of the latter model if one were to neglect 
scattered and re-emitted photons. Similar to the fiducial model 
(Figure [Jjl, scattered and re-emitted photons significantly modify 
the absorption profiles (especially Mg II and Fc II A2600) and pro- 
duce strong emission lines. The result, for Mg II especially, is a set 
of profiles that do not match the average observed LBG absorption 
profile. 

To calculate the optical depth profiles for the other Mg II 
and Fe II transitions, we assume n-pc+ = nyig+ /2 and also 
scale r by the /A product. 

We generated Mg II and Fe II profiles for this LBG- 
Sobolev wind using the ID algorithm with no dust ex- 
tinction; these are shown as red curves in Figure [T5[ For 
this analysis, one should focus on the Mg II A2796 tran- 
sition. The dotted line in the Figure shows the intensity 
profile when one ignores re-emission of absorbed pho- 
tons. By construction, it is the same profile^ described 
by Equation [7] and plotted in Figure 114b . In compari- 
son, the full model (solid, red curve) shows much weaker 
absorption, especially at w = to —300km s~^ because 
scattered photons fill in the absorption profile. In this 
respect, our LBG-Sobolev model is an inaccurate de- 
scription of the observations which show more uniform 
equivalent widths among differing transitions (SIO). The 
model also predicts significant emission in the Mg II lines 
and several of the Fe II* transitions. Emission associ- 
ated with cool gas has beeii observed for Si II* transi- 
tions in LBGs (jPettini et al.ll2002t iShaplev et aIll20Q3| ). 
but the Fe II* transitions modeled here lie in the near- 

^ The other transitions in the LBG-Sobolev model are signifi- 
cantly weaker than the input model because these are scaled down 
by fx and for Fe II the reduced Fe+ abundance. 
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Fig. 16. — The upper left panel shows the density (solid) an d ve - 
locity (dotted) laws for the LBG-Sobolev wind model (Figure 1 1411 . 
The remaining panels show the cumulative mass, energy, and mo- 
mentum of this flow. All of the curves rise very steeply at small 
radii and then rise steadily to large radii. The principal result is 
that the majority of mass, energy, and momentum in the wind are 
contained at large radii. We derive similar results for the LBG- 
partial covering model if we assume the wind is composed of iden- 
tical clouds optically thick to strong metal-line transitions (e.g. 
Si II 1526). The energy and momentum in the outer regions of the 
outflow may be very difficult to generate with standard galactic- 
scale wind scenarios. 

IR and have not yet been investigated. On the other 
hand, there have been no reported detections of signif- 
icant hne-emission related to low-ion resonance transi- 
tions (e. g. Mg II) in LBGs , only z < 1 star-f orming 
galaxies (jWeiiier et all 120091: IRubin et al.1 l2010bD . The 
principal result of the LBG-Sobolev model is that the 
scattering and re-emission of absorbed photons signif- 
icantly alters the predicted absorption profiles for the 
input model. This is, of course, an unavoidable conse- 
quence of an isotropic, dust-free model with unit covering 
fraction. 

Because we have explicit velocity and density profiles 
for the LBG-Sobolev model, it is straightforward to cal- 
culate the radial distributions of mass, energy, and mo- 
mentum of this wind. These are shown in cumulative 
form in Figure [TH Before discussing the results, we of- 
fer two cautionary comments: (i) the conversion of nMg+ 
to nn assumes a very poorly constrained scalar factor of 
10^'^^. One should give minimal weight to the absolute 
values for any of the quantities; (ii) the Sobolev approx- 
imation is not a proper description of the SIO LBG wind 
model (see the following sub-section). These issues aside, 
we may inspect Figure [12] to reveal global properties of 
this LBG-Sobolev wind. One obvious result is that all of 
the curves rise very steeply at small radii. We find, for 
example, that « 10% of the mass is contained within the 
inner 2 kpc. All of the curves continue to rise, however, 
such that the majority of energy, mass, and momentum 
in the wind is transported by its outermost layers (i.e., 
r > 20 kpc) . This is somewhat surprising given that the 
density is « 5 orders of magnitude lower at these radii 
than at r = 1 pc. More importantly, we question whether 
any physical process could produce a wind with such ex- 
treme profiles. 



5.1.2. Partial Covering Fraction 

In the previous sub-section, we described a Sobolev in- 
version that reproduces the average absorption profile of 
LBGs in cool gas transitions when scattered photons are 
ignored. A proper analysis that includes scattered pho- 
tons, however, predicts line-profiles that are qualitatively 
different from the observations because scattered photons 
fill-in absorption and generate significant line-emission 
(similar to the fiducial wind model; §[3]). We also argued 
that the implied mass, energy, and momentum distribu- 
tions of this wind (Figure I16|) were too extreme. This 
LBG-Sobolev model, however, is not precisely the one 
introduced by SIO; those authors proposed an ensemble 
of optically thick clouds with a partial covering fraction 
described by Equation [HI In contrast, the LBG-Sobolev 
model assumes a diffuse medium with a declining density 
profile but a unit covering fraction. One may question, 
therefore, whether these differences lead to the failures 
of the LBG-Sobolev model. 

To more properly model the LBG wind described in 
SIO, we performed the following Monte Carlo calcula- 
tion. First, we propagate a photon from the source un- 
til its velocity relative to line-center resonates with the 
wind (the photon escapes if this never occurs). The pho- 
ton then has a probability P = fc{r) of scattering. If 
it scatters, we track the photon until it comes into reso- 
nance agairFl or escapes the system. In this model, all of 
the resonance transitions are assumed to have identical 
(infinite) optical depth at line center. 

The results of the full calculation (absorption plus 
scattering) are shown as the black curve in Figure 1151 
Remarkably, the results are essentially identical to the 
LBG-Sobolev calculation for the Mg II A2796 transi- 
tion; scattered photons fill-in the absorption profiles at 

V w Okms~^ and yield significant emission lines at 

V > km s^^. Similar results are observed for the other 
resonance lines and significant emission is observed for 
the Fe II* transitions, centered at Sv « 0km s^^. In 
contrast to the Mg II profiles, the Fe II A2586 line 
much more closely resembles its intrinsic profile (dotted 
curves). This occurs because most of these absorbed pho- 
tons fiuoresce into Fe II* emission. The equivalent width 
of Fe II A2586 even exceeds that for Fe II A2600, an 
inversion that, if observed, would strongly support this 
model. 

A robust conclusion of our analysis is that these sim- 
ple models cannot reproduce the observed profiles of res- 
onantly trapped lines (like the Mg II doublet) because 
of the emission from scattered photons. In particular, 
the observations show much greater opacity at veloci- 
ties within « 100km s^^ of systemic (SIO). In order to 
achieve a significant opacity at Sv w Okm s~^, one must 
either suppress the line-emission (e.g. with severe dust 
extinction or anisotropic winds) or include a substantial 
ISM component that absorbs light at 6v > Okm s~^ (e.g. 
§ 14. 3|) . Absent such corrections (which may be insuffi- 
cient), we conclude that the LBG model introduced by 
SIO is not a valid description of the data. Also, similar 
to our LBG-Sobolev calculation (Figure [T4| . this clump 

Because this wind has a monotonically increasing velocity law, 
the photon may only scatter again if a longer wavelength transition 
is available, e.g. Mg II A2803 for a scattered Mg II A2796 photon. 
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model predicts the majority of mass, energy, and mo- 
mentum are carried in the outer regions of the wind. In 
fact, for clumps that are optically thick to strong metal- 
line transitions (e.g. Si II 1526) we estimated!] that the 
wind would carry nearly a total mass of w 10^*^71/0, with 
the majority at large radii. 

The resonant line-emission predicted for Mg II should 
also arise in other transitions. Indeed, we caution that 
the C IV line-emi ssion observed in LBG spectra (e.g. 
iPettini et afl |200^ may result from the galactic-scale 
outflow, not only the stellar winds of massive stars. A 
proper treatment of these radiative transfer effects is re- 
quired to quantitatively analyze these features. 

5.2. Radiation Pressure 

As discussed in the Introduction, the importance of 
various physical mechanisms in driving galactic out- 
flows remains an open question. Several studies ad- 
vocate a primarily energy-driven wind scenario, in 
which the hot wind produced by the thermalized en- 
ergy from SNe eject a entrai ns cool clouds via ram pres- 
sure (e.e. IChevalier fc Clcgd fT985l : fHeckman et all[l99?il : 
[Strickland fc Heckmanii2009l ). Alternatively, momentum 
deposition by photons emitted by luminous star clusters 
onto dust grains in the ISM may also play an integral role 
in dr iving large-scale outflows (e.g. iMurrav et a l. 20Q3, 
l201Clf ). Both of these mechanisms likely contribute more 
or less significantly in a given galaxy; however, the lat- 
ter (radiation-pressure driving) has the added advantage 
that it does not destroy cool gas clouds as it acts. Here we 
explore the line-profiles produced b y a radiation-pressur e 
driven wind model as described bv IMurrav et al.l (|2005[ ). 

These authors start by assuming the galaxy has lumi- 
nosity L and is an isothermal sphere with mass profile 
Mgai{r) = 2a^r/G, where a is the velocity dispersion. 
The gas mass profile is simply Mg{r) — fgMgai{r), with 
fg a constant gas fraction. The dust has opacity k and is 
optically thin, such that the force per unit mass of wind 
material due to radiation is j^^- Ignoring gas pressure, 
and assuming a steady state for the flow (i.e., the mass 
outflow rate, dM^ia^/dt, is constant), the momentum 
equation for the wind is 



-2000 



Velocity (km/s) Relative to Fell 2600 
-1000 1000 



2000 



^dv^^ GMgaijr) 

Substituting 2a^r/G for Mgai{r 
IMurrav et all (pOoBl ) find 



kL 



(10) 



47r cr^ 
and solving for v^, 



v,{r)=2a^\Rg[ - i )+ln(i?oA) (11) 



(their Equ. 26), where v^- — at Rq and Rg — ^^^^'i ■ 
The corresponding density profile for the gas is given by 



dM^ind/dt 



(12) 



with nip the mass of the hydrogen atom. Similar to the 
fiducial wind model, this wind has a decreasing density 

Note that this esti mate differs from the LBG-Sobolev cal- 
culation shown in Figure 1161 because the latter assumes a diffuse 
medium with unit filling factor. 
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Fig. 17. — Mg II and Fe II profiles for a radiation-driven wind 
model (red) compared against the fiducial wind model (black). In 
contrast to the fiducial model, the optical depth peaks at large 
velocity (here &v ft! —350km s"'^). At these velocities, the absorp- 
tion is not fiUed-in by scattered photons and one recovers absorp- 
tion lines that closely resemble the intrinsic optical depth profile 
(dotted curve). Although the Mg II absorption profiles are very 
different from the fiducial model, the Mg II line-emission has sim- 
ilar peak flux and velocity centroids The Fe II* emission also has 
a similar equivalent width, yet the lines are much broader in the 
radiation-driven model. 

law with nn roughly proportional to r~^. In contrast 
to the fiducial model, the velocity law is nearly constant 
with radius before decreasing sharply at large radii. To 
produce a wind whose Mg II A2796 optical depth profile 
peaks at r^gg > 10, we set i?o — ^nnor — 1 kpc, Rg ~ 
4kpc, and a — 125km s~^. The density peaks at r = 
nnncr whcrc — 0.03 cm~'^. The velocity, meanwhile, 
peaks at r = Rg with ~ 400 km s~^. This corresponds 
to a mass flow of dMv^ind/dt w 0.2Mo yr~^. 

Figure [17] presents the Mg II and Fe II profiles for 
this radiation-driven wind model compared against the 
fiducial wind model. In contrast to the fiducial model, 
the optical depth peaks at much larger velocity (here 
Sv w —400 km s^^). At these velocities, the absorp- 
tion is not filled-in by scattered photons and one recov- 
ers absorption lines that more closely follow the intrinsic 
optical depth profile (dotted curve). The Mg II line- 
emission, meanwhile, has similar peak flux and velocity 
centroids to the fiducial model (see also Table [6|). The 
Fe II* emission also has a similar equivalent width, yet 
the lines of the radiation-driven model are much broader. 
This reflects the fact that the majority of absorption oc- 
curs at 6v ^ km s^^. The kinematics of Fe II* line- 
emission, therefore, independently diagnose the intrinsic 
optical depth profile of the fiow. We conclude that the 
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TABLE 5 

Wind Parameters: Power-law Models 



Label 


n(r) 


"h 


v(r) 


vO 


~2796 


A 




0.4000 


^-2.0 


2.0 


0.8 


B 




0.5000 


J.-1.0 


50.0 


1.1 


C 




0.3000 


^0.5 


100.0 


1.8 


D 


^0 


0.0100 


J.-2.0 


2.0 


2.1 


E 


^0 


0.0100 


^-1.0 


50.0 


1.7 


F 


^0 


0.0200 


j.0.5 


100.0 


1.6 


G 


r2 


0.0100 


J.-2.0 


2.0 


4.4 


H 


r2 


0.0010 




50.0 


3.2 


I 


r2 


0.0001 


j,0.5 


100.0 


1.8 


Note. — For velocity laws that decrease with 
radius, the normalization t)" refers to the velocity at 
1" = ''outer instead of r = rjuncr- AH models assume 



a Doppler parameter b = 15 km s ^, full isotropy, 
and a dust-free environment. 

radiation-driven model has characteristics that are qual- 
itatively similar to the fiducial model (e.g. strong blue- 
shifted absorption and significant line-emission) with a 
few quantitative differences that, in principle, could dis- 
tinguish them with high fidelity observations. 

5.3. Power-Law Models 

The fiducial model assumed power-law descriptions for 
both the density and velocity laws (Equations I3l4p . The 
power-law exponents were arbitrarily chosen, i.e. with 
little physical motivation. In this sub-section, we cur- 
sorily explore the results for a series of other power- 
law expressions, also arbitrarily defined. Our intention 
is to illustrate the diversity of Mg II and Fe II pro- 
files that may result from modifications to the density 
and velocity laws. We consider three different density 
laws (nu oc r~'^,r'^,r^) and three different velocity laws 
(ur oc r~^, r""'^, r"'^) for nine wind models (TableEl Fig- 
ure [T8b ) . Each of these winds extends over the same 
inner and outer radii as the fiducial model. The density 
and velocity normalizations n^, have been modified to 
yield a Mg II A2796 optical depth profile that peaks at 
'^2796 ~ 10 — 10'^ and then usually decreases to T2796 < 1 
(Figure ITSb ) . This choice of normalization was observa- 
tionally driven to yield significant absorption lines. All 
models assume a Doppler parameter b — 15 km s~^, full 
isotropy, and a dust-free environment. Lastly, the abun- 
dances of Mg+ and Fe+ scale with hydrogen as in the 
fiducial model (Table [l]). 

Figure [19] presents the Mg II and Fe II profiles for the 
full suite of power- law models (see also Table |6|). Al- 
though these models differ qualitatively from the fiducial 
model in their density and velocity laws, the resultant 
profiles share many of the same characteristics. Each 
shows significant absorption at Sv < km s~^, extending 
to the velocity where T2796 drops below 0.1 (Figure [T8b). 
All of the models also exhibit strong line-emission, pri- 
marily at velocities Sv ~ km s"^. This emission fiUs-in 
the Mg II absorption at 6v > — lOOkms"^ such that 
the profiles rarely achieve a relative flux less than w 0.1 
at these velocities. Similar to the radiation-driven wind 
(i) l5.2p . the power-law models that have significant opac- 
ity at Sv < —300km s^^ tend to have larger peak opti- 
cal depths. Another commonality is the weak or absent 
line-emission at Fe II A2586; one instead notes strong 
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Fig. 18. — In panel (a) we illustrate the density (gray curves; 
dashed - oc r", dotted - oc r~^, dot-dash — oc r^) and 
velocity laws (colored; green - ?; oc r"'^, red - d oc r~^, blue - v oc 
r~^) for a series of power-law models (Table (5)1. For comparison 
we also plot the velocity and density laws for the fiducial model 
(thin, black curves). Panel (b) plots the optical depth profiles of 
the Mg II A2796 line for the nine models built from the density and 
velocity laws of panel (a). Each of these were normalized to have 
peak optical depths of 10 to 1000 and to be optically thin at small 
or large velocity. 

Fe II* A2612 emission that generally exceeds the equiv- 
alent width of the fiducial model. One also notes that 
the width of the Fe II* emission is systematically higher 
for models where the optical depth in absorption peaks 
at large velocity from the systemic. One of the few ob- 
vious distinctions between these models and the fiducial 
wind is the higher peak optical depth of absorption in 
the Fe II A2600 transition. This occurs primarily be- 
cause these models have higher intrinsic optical depths 
(i.e. higher np^+ values). 

While there is commonality between models, the spe- 
cific characteristics of the absorption/emission profiles 
do exhibit significant diversity. Table [6] compares mea- 
sures of the absorption and emission profiles against the 
fiducial model. In detail, the line profiles differ in their 
peak optical depths, the velocity centroids of absorp- 
tion/emission, and their equivalent widths. We find that 
many of the differences are driven by differences in the 
velocity laws, i.e. the kinematics of the outflow. At the 
same time, models with very different density /velocity 
laws can produce rather similar results. For example, 
proflles that have declining opacity with increasing ve- 
locity offset from systemic can be obtained with a radial 
velocity law that increases (e.g. the fiducial model) or 
decreases (model D). The implication is that absorption 
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TABLE 6 

Line Diagnostics for the Fiducial Model and Variants 



Transition 


Model 




Wi 


Wa 




Vt 




^fnt 




/pk 


Vf 


Vf 


Ave 






(km s-l) (A) 






fkm/s") 


("km /si 


Ckm/sl 






(km /s1 


Ckm /s") 


(km /s1 


Mgll 2796 






























Fiducial [- 


-1009.-43 


4.78 


2.83 


0.94 


-215 


-372 


-32,311 


-1.77 


2.48 


32 


117 


215 




LBG Sob. 


-687, -86 


2.62 


1.35 


0.45 


-247 


-347 


-65, 257 


-1.02 


1.84 


21 


90 


193 




LBG Gov. 


-698, -65 


2.67 


1.40 


0.41 


-258 


-344 


-54, 214 


-0.93 


2.23 





72 


172 




Radiation 


-365, -65 


2.86 


1.72 


3.00 


-301 


-250 


-54, 322 


-1.74 


2.13 


53 


119 


236 




A 


-848, -22 


9.29 


5.02 


2.92 


-773 


-499 


-22, 171 


-1.34 


3.30 


21 


64 


118 




B [- 


-1084, -65 


11.08 


7.54 


3.00 


-976 


-589 


-65, 193 


-1.83 


2.89 





50 


172 




G 


-451, -65 


2.54 


1.82 


1.86 


-204 


-200 


-65, 332 


-1.79 


2.09 





116 


247 




D 


[-794, 21 


4.21 


2.14 


1.11 


-11 


-276 


[32,311] 


-1.57 


3.58 


43 


141 


172 


Mgll 2803 




























Fiducial 


-437, -41 


3.29 


1.19 


0.76 


-148 


-193 


-41,676 


-2.24 


2.55 


34 


269 


449 




LBG Sob. 


-426, -73 


1.53 


0.46 


0.24 


-244 


-244 


-62,387 


-0.80 


1.58 


23 


154 


332 




LBG Gov. 


-533, -73 


2.56 


0.98 


0.37 


-233 


-285 


-51,558 


-1.48 


2.35 


2 


227 


396 




Radiation 


-362, -84 


2.32 


1.29 


2.65 


-308 


-256 


-73,323 


-1.31 


1.86 


45 


115 


246 




A 


-554, -41 


7.06 


1.98 


0.78 


-330 


-293 


-41,805 


-5.72 


5.03 


13 


304 


610 




B [- 


-554, -169 


5.74 


1.80 


1.27 


-372 


-359 


-84, 1008 


] -7.31 


4.20 


2 


338 


749 




G 


-405, -73 


2.17 


1.49 


1.85 


-169 


-179 


-73, 355 


-1.52 


1.94 


77 


122 


246 




D 


[-30, 13 


3.69 


0.19 


0.77 


-9 


-1 


[23, 558] 


-1.48 


3.33 


34 


251 


332 


Fell 2586 




























Fiducial 


-348. -35 


0.82 


0.61 


1.01 


-70 


-119 


'-35, 128 


-0.10 


1.11 


35 


47 


128 




LBG Sob. 


-151,-46 


0.18 


0.04 


0.07 


-139 


-99 


[-12.35] 


-0.01 


1.05 





12 


35 




LBG Gov. 


-684, -12 


2.38 


1.93 


0.65 


-162 
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Note. — Listed are the equivalent widths (intrinsic, absorption, and emission), the peak optical depth for the absorption Tpk = — InCJ^min), the 
velocity where the optical depth peaks the optical depth-weighted velocity centroid Vt= = J dv v \n[I (v)] / f dv\n[I(v)], the peak flux /pk in emission, 
the velocity where the flux peaks Vf, the flux- weighted velocity centroid of the emission line Vf (occasionally affected by blends with neighboring 
emission lines), and the 90% width Aue- The v9^^ and vf^^ columns give the velocity range used to calculate the absorption and emission characteristics, 
respecitvely. These were defined by the velcoities where the profile crossed 0.95 in the normalized flux. [See the electronic version for the complete 
Table]. 
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profiles alone are likely insufBcient to fully characterize 
the physical characteristics of the outflow. 

6. DISCUSSION 

We now discuss the principal results of our analysis 
and comment on the observational consequences and im- 
plications. 

The previous sections presented idealized wind mod- 
els for cool gas outflows, and explored the absorp- 
tion/emission profiles of the Mg II AA2796, 2803 doublet 
and the Fe II UVl multiplet. In addition to the ubiqui- 
tous presence of blue-shifted absorption, the wind models 
also predict strong line-emission in both resonance and 
non-resonance transitions. For isotropic and dust-free 
scenarios, this is due to the simple conservation of pho- 
tons: every photon emitted by the source eventually es- 
capes the system to maintain zero total equivalent width. 
A principal result of this paper, therefore, is that galax- 
ies with observed cool gas outflows should also generate 
detectable line-emission. 

Indeed, line-emission from low-ion transitions has been 
reported from star-forming galaxies that exhibit cool gas 
outflows. This includes emission related to Na I (iPhillips 



20101) ■ Mg II emission (jWeiner et al 
2010b( ). Si II* emission from LBGs 



_ _ , Ru bin et al.l . , _ 

IShaplcy et al.l I2003D." and, most r ecently, significant 
Fe II* emission (Rubin et al.|[2010aD . A variety of ori- 
gins have been proposed for this line-emission including 
AGN activity, recombination in H II regions, and back- 
scattering off the galactic-scale wind. Our results sug- 
gest that the majority of the observed line-emission is 
from scattered photons in the cool gas outfiows of these 
star-forming galaxies. Indeed, line-emission should be 
generated by all galaxies driving a cool gas outflow. 

On the other hand, there are many examples of galax- 
ies where blue-shifted absorption is detected yet the 
authors report no s ignificant line emiss i on. This in- 
cludes Na I flo ws (iRu pke et al.l l2005at iMartinI 120051 : 
iSato et al.ll2009D. M g II and Fe II absorption in ULIRGS 
(|Martin fc Bouchel |2009„) . and the e xtrem e Mg II out- 
flows identifled bv iTremonti et ahl (|2007D . Similarly, 
there have been no reports of resonance line-emission 
from low-ion transitions in the LBGfQ, and many z < 1 
galaxies show no detectable Mg II or Fe II* emission de- 
spite significant blue-shifted absorption (Rubin et al., in 
prep). These non-detections appear to contradict a pri- 
mary conclusion of this paper. We are motivated, there- 
fore, to reassess several of the effects that can reduce 
the line-emission and consider whether these prevent its 
detection in many star-forming galaxies. 

Dust is frequently invoked to explain the suppression 
of line-emission for resonantly trapped transitions (e.g. 
Lya). Indeed, a photon that is trapped for many scat- 
terings within a dusty medium will be preferentially ex- 
tincted relative to a non-resonant photon. In § 14.21 we 
examined the effects of dust extinction on Mg II and 
Fe II emission. The general result (Figure [TT|) is a mod- 
est reduction in flux that scales as (1+Tdust)~^ instead 
of exp(-rdust)- Although the Mg II photons are reso- 



IPettini et al.l 112000) reported and analyzed C IV line-emission 
from the lensed LBG cb58, attributing this flux to stellar winds. 
We caution, however, that the galactic-scale wind proposed for 
cb58 should also have contributed to the C IV line-emission. 



nantly trapped, they require only one to a few scatter- 
ings to escape the wind thereby limiting the effects of 
dust. This reflects the moderate opacity of the Mg II 
doublet (e.g. relative to Lya) and also the velocity law 
of the fiducial wind model. In scenarios where the Mg II 
photons are more effectively trapped, dust does suppress 
the emission (e.g. in the ISM-|-wind model from 14.31) . 
In contrast to Mg II, the Fe II resonance photons may 
be converted to non-resonant photons that freely escape 
the wind. Therefore, a wind model that traps Mg II 
photons for many scatterings does not similarly trap 
Fe II AA2586, 2600 photons. The effects of dust are much 
reduced and, by inference, the same holds for any other 
set of transitions that are coupled to a fine-structure level 
(e.g. Si Il/Si II*). In summary, dust does reduce the line 
flux relative to the continuum, but it generally has only 
a modest effect on the predicted Mg II emission and a 
minor effect on Fe II* emission. 

Another factor that may reduce line-emission is wind 
anisotropy. The flux is lower, for example, if one elim- 
inates the backside to the wind (e.g. the source itself 
could shadow a significant fraction of the backside). Sim- 
ilarly, a bi-conical wind can have significantly lower line- 
emission, at least for the fiducial model (§ 14. ip . For 
the anisotropic winds explored in this paper, the line- 
emission scales roughly as fl/Ai: where is the angular 
extent of the wincP^. Because we require that the wind 
points toward us, it is difficult to reduce fl much below 27r 
and, therefore, anisotropy reduces the emission by a fac- 
tor of order unity. Similarly, an anisotropically emitting 
source would modify the predicted line-emission. For ex- 
ample, if the backside of the galaxy were brighter /fainter 
then would one predict brighter/fainter emission relative 
to the observed continuum. To reduce the emission, the 
brightest regions of the galaxy would need to be oriented 
toward Earth. Although this would not be a generic 
orientation, most spectroscopic samples are magnitude- 
limited and biased towards detecting galaxies when ob- 
served at their brightest. In principle, this might imply 
a further reduction of order unity. 

Together, dust and anisotropic models may reduce the 
line-emission of Mg II and Fe II by one to a few factors of 
order unity. For some galaxies, these effects may explain 
the absence of significant line-emission, but they may not 
be sufficient to preclude its detection. 

There is another (more subtle) effect that could greatly 
reduce the observed line-emission: slit-loss. As described 
in Figure |6l emission from the wind is spatially extended 
with a non-zero surface brightness predicted at large 
radii. This implies a non-negligible luminosity emitted 
beyond the angular extent of the galaxy. The major- 
ity of observations of star-forming galaxies to date have 
been taken through a spectroscopic slit designed to cover 
the brightest continuum regions. Slits with angular ex- 
tents of 1 — 1.5" subtend roughly 8 — 15kpc for galaxies 
at z > 0.5. Therefore, a 1" slit covering a galaxy with 
our fiducial wind would cover less than half the wind. 
The result is reduced line-emission when compared to 
the galaxy continuum. 

We explore the effects of slit-loss as follows. We model 

There is also a first order dependence on the actual orientation 
of the wind, with maximal emission when the backside is fully 
viewed. 
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Fig. 19. — Mg 11 and Fe 11 profiles for a series of power-law wind models (see Table[5]for details). These panels reveal the di vers ity of 
absorption and emission that results from a set of intrinsic optical depth profiles that have relatively modest differences (Figure lfSb ). As 
designed, all of the models show significant absorption blueward of systemic velocity. As important, each shows strong line-emission at 
5v ~ km s~^, with the fiux proportional to the degree of absorption. We also note that the width of the Fe II* A2612 emission scales 
with the velocity offset of peak optical depth for the Fe 11 A2586 absorption. Therefore, the emission kinematics trace the intrinsic optical 
depth profile. 



the slit as a perfect, infinitely-long rectangle centered on 
the galaxy. We then tabulate the equivalent width of 
the line-emission through slits with a range of widths, 
parameterized by rrs=o.2'- twice the radius where the 
wind has a Sobolev optical depth of = 0.2. The line- 
emission will be weak beyond this radius because the 
photons have only a low probability of being scattered. 
For the fiducial wind model, the Mg II A2796 transi- 
tion has rrs=o.2 ~ 15kpc (Figure [3]). The predicted 
equivalent width in emission Wo relative to the absorp- 
tion equivalent width M^a is presented in Figure [201 for 
a series of transitions for the fiducial wind model (black 
curves). The We/W^ curves rise very steeply with in- 
creasing slit width and then plateau when the slit width 
reaches ~ ?'rs=o.2- For the fiducial model, the emission 
is concentrated toward the source; this derives from the 
density and velocity profiles but is also a simple conse- 
quence of geometric projection. Figure ?IU\ reveals similar 
results for other wind models. The curves rise so steeply 
that the effects of slit-loss are minor (order unity) un- 
less the slit-width is very small. Nevertheless, the results 



do motivate extended aperture observations, e.g. inte- 
gral field unit (IFU) instrumentation, that would map 
the wind both spatially and spectrally. 

Although several effects can reduce the line-emission 
relative to the absorption of the outflow, our analysis 
indicates that detectable line-emission should occur fre- 
quently. Furthermore, the line-emission could be sup- 
pressed so that it does not exceed the galaxy continuum 
yet still (partially) fills-in the absorption profiles. In- 
deed, dust and anisotropic winds preferentially suppress 
line-emission at 6v > km s~^ fFigures I8I10|) . The re- 
maining emission would still modify the observed absorp- 
tion profiles (e.g. Figure [5]) and may complicate conclu- 
sions regarding characteristics of the flow. An analysis of 
cool gas outflows that entirely ignores line-emission may 
incorrectly conclude that the source is partially covered, 
that the gas has a significantly lower peak optical depth, 
and/or that a 5v ^ 0km s~^ component is absent. We 
now examine several quantitative effects of line-emission. 

Figurc[ni]demonstrates one observational consequence: 
reduced measurements for the absorption equivalent 
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Fig. 20. — Absolute emission equivalent width (|Wc|) relative 
to the observed absorption equivalent width (Wa) for a series of 
transitions: (solid - Mg II A2796; dotted - Mg II A2803; dashed 
- Fe II A2600). The WcjWa. ratio is plotted as a function of slit 
width relative to twice the radius ^,-3=0. 2, defined to be where the 
Sobolev optical depth rg = 0.2. The black curves correspond to 
the fiducial wind model (§|3] ^,-3— 0.2 ~ 15 kpc), the red curves are 
for the LBG-partial covering scenario 15.1.21 r.rs=0.2 ~ 8 kpc), 
and the radiation-driven wind fii 15.21 rVs=o.2 = 40kpc) has blue 
curves. For all of the wind models, the liyol/VVa ratio rises very 
steeply with slit width and then plateaus at fa r^g— 0.2- Therefore, 
a slit that exceeds 0.5rT-g=o.2 will admit nearly all of the photons 
scattered to our sightline. 



Fig. 22. — The observed peak optical depth rpj^ = — ln[/niin] 
for Mg II A2796 profiles of the wind models studied in this paper 
(Tables |4l [G)! with /min the minimum normalized intensity of the 
absorption profile. Cases where exceeds 3 are presented as 
lower limits. The Tpj^ values are plotted against the velocity Vr 
where I(vr) = /min- In of the models, the true peak optical 
depth TpJ,"° > 10. The much lower 'observed' Tp^. values occur 
because scattered Mg II A2796 photons have fiUed-in the absorption 
profiles at velocities 5v > —400 km s~^. These effects, therefore, 
are greatest for wind models where the optical depth peaks near 
5v ~ km because the majority of scattered photons have this 
relative velocity. Indeed, models with Vr < —400 km s^^ all show 
Tpk > 3. 
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Fig. 21. — Difference between the 'observed' absorption equiv- 
alent width Wa which includes the flux of scattered photons and 
the 'intrinsic' equivalent width Wi that ignores photon scatter- 
ing. The dashed (dotted) curves trace a 50% (10%) reduction in 
Wn relative to Wi. One notes a reduction in Wa by 30 — 50% 
for the Mg II A2796 transition (the effect is generally larger for 
Mg II A2803). The effects of scattered photons are reduced for 
the Fe II transitions because a fraction (in fact a majority for 
Fe II A2586) of the absorbed photons fluoresce as Fe II* emission 
at longer wavelengths and do not 'flU-in' the absorption profiles. 
The figure shows results for all of the models presented in Tables |4] 
and|6] 



width Wg, of the flow. In the case of Mg II, which has 
the most strongly affected transitions, is reduced by 
30 — 50% from the intrinsic equivalent width (the equiva- 
lent width one would measure in the absence of scattered 
photons). In turn, one may derive a systematically lower 
optical depth or velocity extent for the wind, and there- 
fore a lower total mass and kinetic energy. The effects 
are most pronounced for wind scenarios where the peak 
optical depth occurs near Sv — km s^^ . Geometric pro- 
jection limits the majority of scattered photon emission 
to have \Sv\ < 200 km s~^; therefore, the absorption pro- 
files arc fillcd-in primarily at these velocities. 

Another (related) consequence is the reduction of the 
peak depth of absorption. In Figure [521 we plot the 
observed peak optical depth Tpk for Mg II A2796 ver- 
sus the velocity where the profile has greatest depth for 
the various wind models (i.e. Tpk vs. Vr from Tables |4] 
and [5]). In every one of the models, the true peak op- 
tical depth Tp™° > 10. For the majority of cases with 

v-r > —300 km s^^, one observes Tpk < 2 and would infer 
the wind is not even optically thick! This occurs because 
photons scattered by the wind have 'filled-in' the absorp- 
tion at velocity Sv « km s^^. In contrast, wind models 
with Vr < —400 km s^^ all yield Tpk > 3. The results are 
similar for the Mg II A2803 profile. In fact, one generally 
measures a similar Tpk for each Mg II line and may incor- 
rectly conclude that the source is partially covered. In 
the case of Fe II, one may even observe a (non-physical) 
inversion in the apparent optical depths of Fe II A2586 
and Fe II A2600 (e.g. Figure [HI). 

These effects are reduced for the Fe II absorption pro- 
files, especially for Fe II A2586. This is because a sig- 
nificant fraction (even a majority) of the Fe II emission 
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FlG. 23. — Comparison of the velocities for which the Fe II A2586 
and Mg II A2796 absorption profiles have greatest depth. Results 
for all of the models studied in this paper are presented (Tables |4l 
[6)l. It is evident that wind models with a peak optical depth near 
the systemic velocity have a Mg II A2796 absorption profile shifted 
blueward by one to several hundred km s^^. Analysis of such 
profiles may lead to the false conclusion that (i) the majority of 
mass in the wind is travelling at a higher velocity; and (ii) there is 
no gas with <5d ~ km s~^ . The dashed curve traces the one-to-one 
line. 

is florescent Fe II* emission at longer wavelengths which 
does not affect the absorption profiles. Therefore, the 
Fe II absorption equivalent widths (Wa) more closely fol- 
low the intrinsic values, one derives more accurate peak 
optical depths, and the absorption kinematics more faith- 
fully reflect the motions of the flow. Regarding the last 
point, one also predicts a velocity offset between the Fe II 
and Mg II absorption- line centroids (Figure [23|). This af- 
fects the analysis of gas related to the ISM of the galaxy 
and also material infalling at modest speeds. Figure [55] 
also emphasizes that the Mg II profiles may misrepresent 
the kinematics of the bulk of the gas. Analysis of these 
lines, without consideration of line-emission, may lead to 
incorrect conclusions on the energetics and mass flux of 
the wind. 

We emphasize that all of these effects are height- 
ened by the relatively low spectral resolution and 
S/N characteristic of the data commonly acquired for 
z > star-forming galaxies. Figure [24k shows one 
reahzation of the Mg II AA2796,2803 doublet for 
the ISM-l-wind model convolved with the line-spread- 
function of the Keck/LRIS spectrometer (a Gaussian 
with FWHM=250km s^^) and an assumed signal-to- 
noise of S/N=7 per 1 A pixel. Both the absorption and 
emission are well detected, but it would be difficult to 
resolve the issues discussed above (e.g. partial covering, 
peak optical depth) with these data. One also notes sev- 
eral systematic effects of the lower spectral resolution, 
e.g. reduced peak flux in the line-emission and a sys- 
tematic shift of Mg II A2796 absorption to more neg- 
ative velocity. We have also modelled an observation 
of the ISM-|-wind model using a stacked galaxy spec- 
trum (Figure [24b). Speciflcally, we averaged 100 iden- 
tical Mg II profiles from the fiducial model degraded to 
a S/N=2pix^^ and shifted by a random velocity offset 
with cr = 100 km s^^. This treatment illustrates the ef- 
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Fig. 24. — 'Observations' of the Mg II profiles for the ISM+wind 
model, (a) The true profile (dotted lines) has been convolved with 
a Gaussian line-spread- function (FWHM = 250 km s~^) and noise 
has been added to give a S/N=7 per pixel. Note the suppressed 
peak in line-emission; one may even infer the flux of Mg II A2796 
exceeds that of Mg II A2803. The Mg II A2796 absorption is also 
shifted to shorter wavelengths, i.e. to a greater velocity offset from 
systemic, (b) The solid curve shows a stack of 100 spectra of the 
ISM+wind profile, each degraded to a S/N=2 per pixel, and offset 
from systemic by a normal deviate with cr = 100 km s~^ to mimic 
uncertainty in the redshift of the galaxies. This treatment is meant 
to illustrate t he implication s of stacking galaxy spectra to study 
outflows (e.g. IWeiner et al.i r2009; Rubin ct al. '2010b, SIO). The 
main difference from the single galaxy observation shown in panel 
(a) is the smearing of line-emission and absorption that reduces 
the height/depth of each. The effects would likely be even more 
pronounced if one studied spectra with a diversity of Mg II profiles. 

fects of using stacked galaxy spectra to stu dy outflows 
(e.g. IWeiner et al. 2009; Rubin ct al.i l2010bl SIO). The 
main difference between this and the single galaxy obser- 
vation shown in panel (a) is the smearing of line-emission 
and absorption that reduces the height/depth of each. 
The effects would be even more pronounced if one stud- 
ied spectra with a diversity of Mg II profiles. Special care 
is required, therefore, to interpret properties of the wind 
(and ISM) from such spectral analysis. 

The previous few paragraphs sounded a cautionary 
perspective on the implications for absorption-line analy- 
sis of galactic-scale outflows in the presence of (expected) 
significant line-emission. While this is a necessary com- 
plication, we emphasize that such analysis remains one 
of the few observational techniques at our disposal to 
study outflows. Furthermore, direct analysis of the line- 
emission offers new and unique constraints on the char- 
acteristics of the outflow. And, when coupled with the 
absorption-line data, the two sets of constraints may 
break various degeneracies in the physical characteristics 
of the outflow. We now consider a few examples. 

The most obvious characteristics probed by the line- 
emission are the si ze and morphology of the outflow 
([Rubin et al.[[2010ah . Line emission is predicted to ex- 
tend to radii where the Sobolev optical depth exceeds a 
few tenths. The principle challenge is to achieve sufficient 
sensitivity to detect the predicted, low surface-brightness 
emission. As Figure [55] demonstrates, the surface bright- 
ness at the inner wind radius of our fiducial wind exceeds 
that at the outer radius by several orders of magnitude. 
Nevertheless, an instrument that sampled the entire wind 
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Fig. 25. — Surface brightness profiles of the fiducial wind model 
for the Mg II A2796 (solid) and Fe II* A2626 (dotted) transitions. 
Both profiles peak at small radii and decrease by several orders 
of magnitude before reaching the outer edge of the wind. The 
profiles arc sufficiently shallow, however, that the azimuthally in- 
tegrated flux declines by only a factor of as 10 from r = 1 kpc to 
r = 10 kpc. Therefore, a sensitive IFU observation could map the 
emission (and, in principle, the kinematics) from rinnor to approx- 
imately router- 

(e.g. a large format IFU or narrow band imager) may de- 
tect the emission in azimuthally-averaged apertures. For 
the fiducial wind model, for example, the azimuthally in- 
tegrated flux falls by only a factor of 10 for projected radii 
of 1 to 15 kpc. Of the transitions considered in this pa- 
per, Mg II emission is preferred for this analysis because 

(i) it has the largest equivalent width in absorption and 

(ii) there are fewer emission channels per absorption line 
than for Fe II. The Mg II lines (especially Mg II A2803) 
frequently have the highest peak and integrated fluxes 
(Figure [25]) . On the other hand, the Mg II transitions 
are more susceptible to the effects of dust extinction and 
for some galaxies Fe II* emission could be dominant. In 
either case, the study of spatially-extended line-emission 
from photons scattered by a galactic-scale wind offers a 
direct means to study the morphology and radial extent 
of these phenomena. The next generation of large- format 
optical and infrared integral-field-units are well-suited to 
this scientific endeavor (e.g. KCWI on Keck, MUSE on 
the VLT). 

The kinematic measurements of the line-emission also 
offer insight into physical characteristics of the wind. In 
§ 14.21 we emphasized that dust extinction preferentially 
suppresses photons scattered off the backside of the wind 
(e.g. Sv > km s^^) so that the line-centroid is shifted to 
negative velocities. Similarly, anisotropic models with re- 
duced emission off the backside yield emission lines that 
are centered blueward of the galaxy's systemic redshift. 
These effects are most prominent in the Fe II* emission. 
The centroids of the Mg II lines, meanwhile, are sensi- 
tive to the optical depth profile of the wind. For ex- 
ample, a wind with fiows exceeding « 600km will 
shift the Mg II A2796 centroid to bluer wavelengths (e.g. 
Figure fT3|) . Although these effects are modest (tens of 
km s~^), they may be resolved by moderate resolution 
spectroscopy. 

The emission-line velocity widths also reveal charac- 
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Fig. 26.— Plot of the 90% velocity width (Av^) of Fe II* A2612 
emission versus the velocity where the resonance lines achieve peak 
optical depth Vt in absorption. The Ave values rise steadily with 
increasing offset of Vt from systemic. The width of the Fe II* A2612 
emission line, therefore, offers an independent diagnostic of the 
wind speed. 

teristics of the wind. In particular, the Fe II* emission- 
line widths are sensitive to the optical depth of the wind 
at large velocity offsets from systemic. This is illus- 
trated in Figure [221 which plots the 90% velocity width 
Ave of Fe II* A2612 versus the velocity where the wind 
optical depth is maximal (wr)- We find that Ave in- 
creases with Vr such that a large width for Fe II* emis- 
sion requires a wind profile with large optical depths at 
Sv < —200 km s^^. The broadening of Fe II* emission 
occurs because this emission is dominated by single scat- 
terings which trace all components of the wind. 

Lastly, the line flux ratios of pairs of transitions are sen- 
sitive to characteristics of the wind. The most obvious 
example is the relative emission of the Mg II AA2796, 2803 
doublet. Specifically, an outflow whose velocity exceeds 
the doublet spacing {Av « 770km s^^) will convert 
Mg II A2796 photons into Mg II A2803 emission resulting 
in a flux ratio that is inverted relative to the intrinsic op- 
tical depth proflles (e.g. Figure [T3l) . Another important 
example is the relative flux of Fe II* emission relative to 
Mg II. A dusty medium with velocity near systemic (e.g. 
an ISM component) may signiflcantly suppress Mg II 
emission yet have a significant Fe II* line fiux because 
the latter is not resonantly trapped. 

7. SUMMARY 

In this paper, we have explored the predicted absorp- 
tion and emission-line profiles of a set of simple galactic- 
scale outflow models. This analysis implemented Monte 
Carlo radiative transfer techniques to propogate reso- 
nant photons through an expanding medium, allowing 
for their conversion to non-resonant photons (e.g. Fe II*). 
Our work focused on the Mg II AA2796, 2803 and Fe II 
UVl multiplet of rest-frame UV transitions, but the re- 
sults apply to most other lines used to probe cool gas. 
Our primary findings are summarized as follows: 

1. Isotropic, dust- free wind models conserve photon 
flux. Therefore, the blue-shifted absorption-line 
profiles commonly observed in star-forming galax- 
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ies are predicted to be accompanied by emission- 
lines with similar equivalent width. This holds even 
for non-extreme anisotropic and dust-extincted sce- 
narios. 

2. The line-emission occurs preferentially at the sys- 
temic velocity of the galaxy and 'fills-in' the absorp- 
tion profiles at velocity offsets \6v\ < 200km s~^ 
from systemic. For transitions that are only cou- 
pled to the ground-state (e.g. Mg 11, Lya, Na I), 
this implies much lower, absorption-line equivalent 
widths (by up to 50%) and observed absorption 
profiles that are significantly offset in velocity from 
the intrinsic optical depth profile. 

3. Analysis of cool gas outflows that entirely ignores 
this line-emission may incorrectly conclude that the 
source is partially covered, that the gas has a signif- 
icantly lower peak optical depth, and/or that gas 
with velocities near systemic (e.g. from the ISM or 
even an infalling component) is absent. 

4. Resonance transitions that are strongly coupled to 
non-resonant lines (e.g. Fe II, Si II) produce emis- 
sion dominated by the optically-thin, fine-structure 
transitions. As such, these resonance absorption 
lines offer the best characterization of the opacity 
of the wind and also of gas with velocities near sys- 
temic. 

5. Dust extinction modestly affects models where res- 
onance photons are trapped for only a few scatter- 
ings. Models with high opacity at small radii and at 
systemic velocity (e.g. with an optically thick, ISM 
component) can effectively extinguish resonantly 



trapped emission (Mg II) but have weaker effect 
on non- resonant lines (Fe II*). 

6. We examined two scenari os designed to mimic the 
wind model proposed by 'Stei del et al.l (|2010[ ) for 
z ~ 3 Lyman break galaxies. Our implementation 
of this model genreates substantial line-emission 
from scattered photons that greatly modifies the 
predicted line-profiles so that this model does not 
reproduce the observed line-profiles. 

7. Significant line-emission is a generic prediction of 
simple wind models, even in the presence of dust, 
anisotropic flows, and when viewed through finite 
apertures. We have explored the 2D emission maps 
(Figure |6]) and surface brightness profiles (Fig- 
ure [25]) of the winds. Sensitive, spatially-extended 
observations will map the morphology and radial 
extent of the outflows. These data afford the best 
opportunity to estimate the energetics and mass- 
fiux of galactic-scale outfiows. 

8. The kinematics and fiux ratios of the emission 
lines constrain the speed, opacity, dust extinction, 
and morphology of the wind. When combined 
with absorption-line analysis, one may develop yet 
tighter constraints on these characteristics. 
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